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Abstract 


A  previous  study  at  the  Air  Force  Institute  of 
Technology  has  led  to  the  development  of  a  closed  form 
solution  to  the  dynamic  behavior  of  a  charged  particle  beam. 
This  solution  was  found  for  a  simple,  but  nontrivial,  single 
degree  of  freedom  model  (the  "electrostatic  approximation 
model").  This  study  investigates  the  characteristics  of 
this  closed  form  solution  using  three  different  types  of 
inputs.  The  closed  form  solution  is  then  verified  by  using 
finite  differences  to  generate  compatible  results.  Also 
investigated  are  the  characteristics  of  the  dynamic  model 
prior  to  the  implementation  of  the  electrostatic 
approx imat ion. 


xi 


CHAPTER  1 


1 -  1  Introduct Ion. 

In  Reference  [1]  a  system  of  partial  differential 
equations  is  presented  as  a  description  of  the  dynamic 
behavior  of  a  particle  beam.  This  system  is  then  cast  as  an 
abstract  Cauchy  problem.  This  results  in  a  set  of  equations 
in  the  following  form: 

<d/dt)  X  =  AX  +  g(t)  (1.1) 

Where  A  is  a  non-constant  9x9  matrix  and  cj(t)  is  the 
input  vector. 

In  til,  this  problem  is  analyzed  using  semigroup  theory, 
from  which  a  closed  form  solution  is  derived  for  a  single 
degree  of  freedom  model  (the  "Electrostatic  Approximation 
Model") . 

In  this  thesis,  numerical  solutions  are  obtained  using 
this  closed  form  solution.  Also  an  alternative  solution  of 
the  Electrostatic  Approximation  Model  is  derived  as  is  an 
approximate  solution  of  the  full  system.  See  Figure  1  for  a 
block  diagram  representation  of  the  overall  flow  of  this 
thesis.  Section  1-2  describe  in  more  detail  the  contents  of 
the  following  chapters. 


9x9 

Linear 

Model 

x  ( t ) 

=  A  x  ( t ) 

+  g<  t  > 

T 


1-2  Chapter  Descriptions 


In  Chapter  2  of  this  thesis  the  closed  fora  solution 
discussed  in  Section  1-1  is  investigated  using  three 
distinct  types  of  inputs.  The  first  input  is  a  step 
function  which  is  Implemented  at  t  =  0.0.  The  second  input 
is  a  function  which  exponentially  approaches  its  maximum 
value.  The  third  input  rapidly  reaches  its  maximum  value, 
is  constant  for  an  interval  of  time  and  then  rapidly  returns 
to  equ 1 1 Ibr l urn . 

In  Chapter  3,  the  system  of  equations  discussed  in 
Section  1-1,  Equation  (1.1),  are  solved  numerically  after 
the  system  is  reduced  to  fifth  order  by  means  of  the 
Electrostatic  Approximation  Model.  By  approximating  the 
partial  derivitives  in  r  with  finite  difference  quotients 
and  then  solving  the  resulting  differential  equations,  the 
solutions  are  determined.  Since  four  of  the  five  equations 
in  this  system  do  not  involve  derivatives  with  respect  to  r, 
they  are  solved  exactly.  The  input  used  in  this  chapter  is 
the  step  input  described  earlier. 

The  analysis  presented  in  Chapter  4  is  closely  related 
to  the  Chapter  3  material  with  the  Electrostatic 
Approximation  Model  not  being  used.  Therefore,  the  system 
shown  in  Section  1-1,  Equation  (1.1)  is  of  ninth-order. 

Also,  this  system  of  differential  equations  are  even  more 
complex  since  five  of  the  nine  equations  involve  finite 


difference  approximations.  The  same  Input  to  the  system 
used  in  Chapter  3  is  also  used  for  Chapter  4,  however  the 
solution  technique  is  considerably  different  due  to  the 
complexity  of  the  problem.  This  ninth  order  system  of 
differential  equations  is  solved  by  an  IMSL  routine  called 
DGEAR  using  the  Cyber  Computer  System. 

Chapter  5  of  this  thesis  contains  conclusions  and 
recommendations  along  with  all  comparisons  of  data  and 
discussion  of  results. 


CHAPTER  2 


2-1  Introduct Ion. 

The  purpose  of  this  section  is  to  use  the  closed  fore 
solution  to  the  reduced  order,  linear,  single  degree  of 
freedoa  Model  found  in  Reference  il:  page  IV-211  for  a  given 
input  and  calculate  the  perturbed  state  of  the  particle  beaa 
in  question.  The  results  fro*  this  section  are  used  later 
for  coaparison  with  results  froa  other  procedures. 

The  particle  beaa  used  throughout  this  effort  is 
assuaed  to  be  at  equillbriua  for  t  <  0.  In  order  to  aeet 
this  assuaption,  an  external  Magnetic  field  <B0>  with  only 
the  Z  coaponent  non  zero  is  applied  for  all  tlae  (t).  The 
calculations  required  to  deteraine  B0  are  found  in  Section 
2-3.  The  beaa  paraaeters  are  as  follows: 

Current  (I)  ■  1  aaperes  (asps) 

Radius  (Rb>  -  10  allliaeters  (aa) 

Kinetic  Energy  CKE)  *  10  alllion  electron  volts  (aev) 
Hass  <M0)  *  1.6726825  x  10'27  kilograas  (kg) 

Charge  (q)  *  1.60219  x  10" 19  couloab 


In  Reference  (1:  page  1—11  a  set  of  partial 


differential  equations  that  describe  the  dynaaic  behavior  of 
a  particle  beaa  are  cast  into  the  fora  of  an  abstract  Cauchy 
problea.  The  ninth-order  linear  systea  of  equations  is  of 
the  fora 


<  d/dt  >  =  A  g<t>  +  g< t )  (2.  la) 

where  £(t)  is  the  perturded  values  froa  equilibrium  and 
g(t)  and  A  are  defined  as  follows: 
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*Ere<r,t>  -  W  r  $BzeCr, t)  -Vz°  *Bee(r,t> 
dEge (r , t )  +  Vz°  #Br«Cr,t> 

*Ez®<r,t>  +  «  r  #Br*<r,t> 


( 2 . lb) 


0 


0 


**25  vz°  “*25  w  r 


B  i 

■  r 


*25 

0 

0 

0 

0 

0 

-Dr 


0  0  0 

«25  0  0 

0  0  0 

0  0  -C2  Dr 

0  C2 

Dr  0  0 
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where 


*61  "  *  <  *52  W  r  )  /  M° 
*71  "  <  *52  vz°  >  f  M° 


(2. 1c) 


o  * 


and  a23»  a24»  a25»  »52»  Dr»  Dr»  are  referenced  In  the 
syabols  page  and  defined  In  Equations  (2.7a)  to  (2.7k) 
Next,  this  system  of  equations  can  be  reduced  by  the 
implementation  of  the  Electrostatic  Approximation  Model 
described  in  Chapter  1  and  fully  developed  in 
Reference  (1:  page  1 1 1-391.  This  approximation  reduces  A 
and  g(t)  of  Equation  (2.1)  to  the  following  much  more 
manageable  system. 


A  * 


0 

0 

0 

0 

0 


-N°  Dr 
0 
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~®24  r 
a52 


0 
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a24  r 
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(2.2a) 


g(t) 


0 


a25 


dEr* (r , t)-W  r  *Bza(r,t>-Vz°  *Be«(r,t) 
#Ee*(r,t)+Vz°  *Br*(r,t> 

$Bze (r, t)+W  r  #Br*(r,t> 


0 


(2.2b) 


The  closed  fora  solution  found  in  Reference  (1:  page 


IV-211  to  this  fifth  order  system  of  equations  and  aentloned 
in  the  previous  section  is  as  follows. 


X<t)  *  S(t)g°  +  \sct-s)g(s>ds  (0  <  t  <  T) 


where 


(2.3) 


}(<t>  » 


«N(r,t> 
#Vr<r,t> 
dVeCr.t) 
#V2<r,t> 
*Er(r, t) 


g(s)  is  as  defined  in  Equation  <2. 2)  where  t  Is  replaced 
by  s*  and,  from  Reference  11:  page  IV-221,  the  state 
transition  matrix  S<t>  Is  the  following  5x5  system 
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M°  Dr  slntQtK.) 


M°  Dr  (cos  (Qtl-1 )  «23  <•> 
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(2.4) 


(costQt l-l ) 


*25  *52 

l  -  (costQt 1- l ) 


q2 


-10- 


. *»' 4 % V- * 


£°  Is  the  initial  state  of  X(t>  and  hence  is  the  difference 
between  the  initial  state  of  the  beaa  and  the  equil ibrlua. 
Since#  in  the  problem  stated  above,  these  quantities  are 
equal , 

X°  ■  0  (2.5) 

Equation  (2.3)  then  becoaes 

( 

25< t)  -  \  S(t-s)g(s)ds  for  (0  <t  <T)  (2.6) 

O 

It  should  be  noted  that  the  external  magnetic  field 
required  to  achieve  equilibrium  is  not  included  in  the  input 
tern  g(a).  This  is  because  the  formulation  of  the  problem 
leads  to  the  solution  of  Equation  (2-3)  in  which  perturbed 
quantities  depend  only  on  perturbed  inputs  #Er*  through 

fBj.®  in  g(s)  while  equilibrium  quantities  appear  in  S(t). 


2-3  Calculation  of  the  Entries  of  S(t-s)  and 


Now,  to  calculate  the  entries  of  S<t-s)  and  g(s)  of 
Equation  (2.6),  t,  Vz°*  N0,  ^p*  a24*  a25»  a52  ®o»  ^c»  w* 
a23»  and  Q  must  first  be  determined.  These  quantities  are 
defined  in  the  list  of  symbols  pages  viii  to  x  and  are  now 
numerically  calculated. 

2.3.1  Calculations  of  t,  V2°,  N°,  Wp 

From  Reference  (2:  page  21,  if 

K.E. 

t  =  1  +  -  (2.7a) 

R.M.E. 

=  1.01065 


then 

V2°  =  (  1  -  1/t2  )1/2  C  (2.7b) 

=  4.34198  x  107  */sec 

For  the  remainder  of  this  thesis  the  beam  is  assumed  to  be 
non-re lat l v ist 1c  and  t  Is  assumed  to  be  equal  to  one. 


From  Reference  £2:  page  11  assuming  uniform  current  density 


-  4.57559  x  1014  l/»3 
Fro*  Reference  (2:  page  3] 


W  2 


<  UQ  C2  M°  ) 


so 


Wp  ■  2.81613  x  107  rad/sec 


2-3.2  Calculations  of  ag4»  ags*  *52*  Bo 


Pros  Reference  11:  page  I I 1-371 


*24 


-Wp2  Vz° 


2  C* 


-1.91571  x  105  <a-sec>"!  , 


*25  “ 


9.57856  x  107  couloab/kg 


*52  *  “°o  9  c2  1,0 


(2. 7g) 


■  -8.27953  x  to6  kg/sec2-couloab 
Froa  Reference  (2:  page  3) 


2  M0  I 

q  t  if,  V2°  Rb2 


<2.7h) 


so 


B0  *  .41578  tesla 


2-3.3  Calculations  of  Wc(  W.  *23*  Q 


Froa  Reference  (2:  page  31,  Wc  is  calculated  by 


q  B0 

Wc  ■  -  (2.71) 

Mo 


■  3.98261  x  107  rad/sec  , 


Froa  Reference  (2:  page  31 


«0  ■  -Wc  <1  +  <1  -  2  Wp2  /  Wc2)l/2>  /2 


<2 ,7j 


It  is  iaportant  to  note  here  that  W$  found  above  is 


related  to  W  of  Reference  tl:  page  I I 1-27 1  by  the 
following  relationship. 

W  -  -W0 

The  value  for  W  is  used  froa  this  point  on.  This  is 
done  to  comply  with  the  nomenclature  found  In  Reference 
tl:  page  I I 1-371. 

In  Reference  tl:  page  111-37],  a23  i*  defined  to  be 

*23  ■  Wc  -  2  W  <2.710 

One  type  of  flow  which  is  often  considered  is  Brtllouln  flow 
<s«e  Reference  (2:  page  4:  case  2J  >.  For  Brillouln  flow 
and  using  the  definitions  found  above,  H  is  calculated  to 
be 

W  -  »c  /  2  <2. 71) 

*  19913067.12803  rad/sec 

Thus,  for  Brlllouin  flow 

023  *  0*0  rad/sec 

A  key  frequency  which  is  used  later  is  denoted  by  the  symbol 


fl  and  is  defined  by 


Q  =  C  Wc2  -  Wp2<  1  -  A  )  i1/2  <2. 8) 

where  :  A  =  Wp2 

For  the  values  of  r,  Vz°,  and  C, 

A  is  sufficiently  close  to  zero  to  be  considered 
negl igible . 

Moreover,  Wc2  =  2  Wp2. 

Therefore  Equation  (2.8)  becomes 

Q  %  Wp 

=  2.81613  x  107  rad/sec  2-4 


i(s)  .S(t-s)< 


and  X(t) 


At  this  point  all  the  terms  necessary  to  completely 
determine  a  solution  to  Equation  (2.1)  starting  with  A  and 
g<s)  defined  in  Equation  (2.2)  are  available. 

To  maintain  the  rigid  rotor  equilibrium,  the  only 
external  field  which  is  non-zero  is  the  Z  component  of  the 
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external  Magnetic  field  <B0>.  For  this  thesis,  the 
perturbed  Input  Is  a  small  perturbation  of  this  external 
field  and  consequently  is  taken  to  be  one  percent  of  B0. 
Therefore  the  input  g(s>  reduces  to  : 


g<s>  -  a23 


0 

-W  r  *Bz*Cr,s> 

0  ( 
0 
0 


*  *23 


0 

C  -82795.32170  kg/couloab  -  sec2  )  r 
0 
0 
0 


The  next  step  is  to  find  the  product  of  8(t-»>g(s>  and 
perfora  the  required  Integration. 


Using  ths  definition  of  8<t-s>  and  g(s>: 


S(t-s)g(s>  -  a2S 


~H°  Dr  ainl  )  r  <  D  ) 

Q 

cost  Q(t-s)  1  r  <  D  ) 


-023  r  slot  QCt-s)  )  (  D  > 


-o24  r 


2  sint  Q( t-s>  1  (  D  > 


052  r  sint  OCt-s>  ]  (  D  ) 


(2.10) 


where  0  »  -82795.32170  kg/coulonb-sec2 


Nos,  to  perfors  ths  rsquirsd  integration  necessary  to 
deterslne  JJ(t)  in  Equation  (2.7),  the  tr igonosetrlc 
functions  of  Equation  (2.10)  can  be  expanded  by  angle  sua 
and  angle  difference  relations.  After  thts  expansion  is 
completed  and  constants  are  factored  out  of  the  integral  the 
resulting  values  can  be  determined.  They  are  as  folloms: 


-2  D 


*<t>  -  *25 


O2 
D  r 


(  1  -  cost  fit  )  > 


-  D  *23  r 

q2 

-  D  *24  r2 

fi2 

~  P  *52  r 
02 


sint  Qt  1 


<  1  -  cost  fit  J  > 


<  1  -  cost  fit  I  > 


<  1  -  cost  fit  )  > 


<2.11) 


and,  substituting  the  appropriate  values  for  the  constants 
gives: 


X<t> 


<  9151.20373  1/aa3  )  *  <  1  -  cost  fit  1  > 

<  -281613.48000  1/sec  >  *  r  *  sint  fit  1 

<  0.00  > 

<-1915.71291  1/a-sec  >  *  <  r2)>  *  <  1  -  cost  Qt  1> 

<  82795.44549  kg/sec2-couloab  )*<  r  >*<  l-cost  fit  )> 


The  above  results  are  tabulated  In  Tables  1  thru  4  and 


are  shown  graphically  in  Figures  2  thru  5 


In  this  section  the  work  presented  In  Section  2-4  Is 
repeated,  only  this  tiae  using  different  external  inputs 
dB2e(r,t).  The  first  of  these  inputs  is  designed  to  aore 
realistically  aodel  how  an  actual  Input  would  be  applied  to 
the  systea.  In  other  words,  the  input  reaches  its  aaxlaua 
value  rapidly  instead  of  instantaneously  as  assuaed  in 
Section  2-4,  and  as  Indicated  in  Figure  18.  This  work  is 
presented  In  Section  2-5.1.  The  second  input  explored  In 
this  section  is  a  fluctuation  froa  the  equiltbrlua  value  of 
*V»<r  ,t>  which  rapidly  reaches  its  aaxlaua.  Then,  after  a 
specified  aaount  of  tiae  <t|)  it  decreases  back  to  its 
equillbrlua  value  (see  Figure  19  >.  The  input  for  this  third 
case  is  the  saae  as  that  in  Figure  18  for  0  1  t  i  6  x  10“09 
seconds.  It  Is  then  constant  for  6  x  I0“M  1  t  1  444  x 
10“09  seconds.  The  Input  then  decreases  as  indicated  In 
Figure  19  for  444  x  I0"0S  <  t  <  450  x  10'°®  seconds.  For  t 
>  450  x  10~09  seconds  the  input  is  zero. 


2-5. 1  BflBlfl_bai-£2QUaS2aB_IaSESftB*_i2_£2QfiSfiBl_ 

Using  the  new  input  described  In  Section  2-5  the  g(s) 
•atrtx  becoses: 

0.0 

-W  r  #BX® 

g(s)  *  *25  0.0 

0.0 

0.0 

where 

#Bz«<r,s>  -  .01  B°  (  t  -exp(  -H  s  )  ) 

and 

H  «  1  x  1009 

The  dynaalc  behavior  of  this  new  Input  Is  shown  In  Figure 
18. 

The  last  step  required  before  Integration  Is  to  find 
the  product  of  the  3<t-s)  and  g<s)  aatrloes,  where  the 
S<t-s>  uatrlx  Is  exactly  the  sane  as  used  In  Section  2-4  and 
the  g<s>  aatrlx  Is  as  defined  above.  This  resulting  product 


9(t-s)g(s>  -  «25 


N<>  Dr  slnt  Q<t-s)  1  C  r  W  #BZ® <  r ,  s  >  ) 
Q 

-cost  Q(t-s)  1  W  r  #Bze<r,s> 

*23  stnt  Q(t-s)  i  W  r  #Bz*<r,s> 


Q 

*24  r  slnt  Q(t-s)  i 

W  r  #Bze<r,s> 

0 

-*52  sint  Q(t-s)  1 

«  r  #Bz*<r,s> 

Now,  rearranging  the  terns  of  the  S(t-s)  g<s)  matrix  yields 


S(t-s>  g(s>  ■  *25 


2  W  NO  tint  0<t-s>  1  #Bz*<r,s> 
fi 

-  W  r  cost  0Ct-s>  1  #8z*(r,s> 

*23  «  r  Btnl-fllt-I?.  )  «z*Cr,s) 
0 

*24  r2  w  alnl  .fln-g?,  l  «z*cr,s> 
0 

-•52  «  r  slnt  0<t-s)  1  #Bz*<r,s> 
P 


Before  perforelng  the  required  integration  on 
8<t-s>g<s>  it  is  important  to  notice  that  this  requires 
the  integration  of  a  product  of  the  fore: 


slnt  P(t-s)  1  ezp<  -a<s>  > 


This  type  of  integration  can  be  performed  by  first 


expanding  the  trigonetric  functions  using  angle  sub  and 
angle  difference  relations  and  then  integrating  the 
individual  terns  using  standard  integration  foraulae. 
Therefore,  the  results  for  Equation  (2.6)  are  as 
follows: 


JtKr,t> 


<  .02)  B®  a25  W  N°  (  X  +  Y  -  Z  ) 
Q 


(2.12) 


where 


X  »  1  -  cosflt 

0 


exp<  -H  t  ) 


H  slni  2 lit  1  ♦  Q  cost  20t  j] 
L  <  ~H  >2  +q2 J 


Z  - 


JVr(r.t)  -  -a2g  W  r  (.01)  B°  (  X2  ♦  Y2  ♦  Z2  >  (2.13) 


*h*re  X2  - 

Q 


Y2  - 

Z2  - 


tLilPi.-H 

(  -H)'z  +  Q* 


± _ CQS  [  Qt  1  -ft  sin  (  Qt  1 

(-H)"  +  Q2 


#V#(r,t>  -  (.01)  B°  a23  a2g  W  r  (  X3  ♦  Y3  -  Z3  >  (2.14) 


where  X3 


Q 


<  -H  >z  +  Qz 


a24  >25  r2  w  <-01)  B°  C  X4  +  Y4  -  Z4  )  (2.15> 

#Vz(r,t>  -  - 

O 


z<  -  i-±.fl-sg>L.e$.  i- 

<  -H  >z  +  Q* 


#Er(r>  t >  -  -a52  *25  «  r  C.01)  B°  (  X5  +  Y5  -  Z5  )  (2.16> 

O 


The  results  from  this  section  are  tabulated  in  Tables 
5  thru  I  and  are  shown  graph  leal If  in  Figures  6  thru  9. 


The  third  and  final  type  of  external  input  £Bze(r,t) 


explored  in  this  chapter  is  a  fluctuation  away  from  zero 
which  returns  to  zero.  The  procedure  presented  in  8ectlon 
2-5.1  is  repeated  using  this  new  input.  Since  the 
procedures  are  nearly  Identical*  the  work  presented  in  this 
section  is  without  explanation.  If  a  question  arises  from 
the  analysis*  refer  to  Section  2-5.1  for  an  explanation. 


The  input  Is  defined  by 

0.01  B°  <l-e"H  *> 


0<t  <t|  <2. 17a) 


«B,#<r.t>-  ^ 0.01B° 

V  »•« 


1-  3<t-ti>2  ♦  2Ct-tt>3 


<t2-t|)2  (tj-tj)3 


tt<t<t2  <2. 17b) 

t)t2 


where  tj«4.45E-7 
t2  »  4.50  K  -  7 

Using  this  definition  for  the  external  input  requires 
that  the  Integration  of  the  S(t-s)  g<s)  matrix  be  perforaed 
over  two  Intervals  of  tine  with  each  Interval  having  a 
different  Integrand  which  involves  the  external  input  as 
shown  in  the  above  equation.  The  Integration  over  the  first 


interval  has  already  been  done  in  Section  2-5.1.  Therefore , 
all  that  remains  is  to  perform  the  integration  over  the 
second  tine  interval  and  add  the  solution  found  in  Section 


2-5. 1.  The  integration  of  the  S(t-s>  g<s>  matrix  over  the 
second  time  Interval  results  in  the  following: 


JN<r,t>  ■  0.02  B0  a2s  «  H°  i  -  cositHt-tj)) 


(2. 18) 


<t2-tj )2 


Q2  <t-t|>  -  2  +  2  coslQ(t-tj)) 


ftt2  <t-tt>3  -  6  Ct-tj)  ♦  6  s ln(fi< t-t | > I 


dVr(r,t>  »  -0.01  B0  a2g  W  r  sin(  fc<t~t|>  ) 


C2.19) 


3  2  <t-tt>  -  2  sinCQ(t-t |) i 

< t2-t i )2  q2  q3 

2  [3  a2  <t-t|)2  -  6  ♦  6  coslO<t-t|>i 

<t2-tj>3  Q4  Q4 


«V#<r,t>  -  0.01  Bq  a23  a2g  If  r  ["l  ■  cosCO(t-tj))  (2.20) 


2 


Over  the  third  interval.  t2  <  a  <  t,  the  Input  la  defined 
by  Equation  (2.17>  aa  zero  and  the  aolutlon  la  the  auaaation 
of  the  Section  2-9.1  aaterlal  Integrated  from  t  ■  0.0  to  t  ■ 
t{  with  Squat  Iona  <2.18)  thru  <2. 22)  integrated  from  t  ■  tj 
to  t  -  t2.  The  reaulta  are  preaented  In  Tables  9  thru  16 
and  are  shown  graphically  In  Figures  10  thru  19.  They 
represent  the  results  froa  this  section  after  each 
particular  type  of  input  was  applied  to  the  systea  for  Its 
associated  aaount  of  tlae. 


CHAPTER  3 


3- 1  Introduction. 

This  chapter  presents  an  analysis  of  the  linear,  time 
invariant,  partial  differential  equations  that  govern  the 
dynamic  behavior  of  the  particle  beam  discussed  in  Chapter 
2.  However,  In  this  chapter  the  solution  to  the  governing 
equations  is  found  numerically.  This  Is  done  by 
approximating  the  partial  dertvltives  In  r  with  finite 
difference  quotients  and  then  evaluating  the  equations  along 
the  radius  of  the  beam  at  selected  points.  These  results 
are  tabulated  and  compared  to  the  results  found  in  Chapter  2 
from  the  closed  form  solution  technique  derived  in  Reference 
Cl:  page  I V—  211. 

3-2  Setting  Up  the  Problem. 

From  Reference  Cl:  page  1 1 1-401  the  governing  partial 
differential  equations  discussed  above  are  as  follows  : 


<  d/dt  )  X(t)  = 

0  -N°  Dr  0  0  0  " 

0  0  a23  a24  r  a25 

0  -a23  000  X(t>  (3.1) 

0  -a24  r  0  0  0 

.0  a52  0  0  0  _ 


0 

*Er«(r,t)  -  W  r  tfBze(r,t>  -  Vz°  *Be«(r,t 
+  a25  $Eee(r,t)  +  Vz°  $Bre(r,t> 

£Eze  <r ,  t)  +  W  r  SBre (r ,  t) 

0 


where  Dr  f  =  <l/r)  Dr<  r  f(r)>  (3.2) 

and  where  the  values  for  a23,  a24,  a25>  and  ag2  are 
the  same  as  defined  in  Chapter  2.  Also,  as  in  Chapter  2, 
the  only  external  field  applied  is  £Bze(r,t)  and  its  maxiBum 
value  is  assumed  to  be  one  percent  of  the  equilibrium  value 
B0 .  The  rest  of  the  external  fields  are  assumed  equal  to 
zero.  Therefore,  Equation  (3.1)  becomes  : 


<  d/dt  )  X(t)  = 


0  -N°  Dr  0 

0  0  a23 

0  -a23  0 

0  “a24  r  0 

0  a52  0 

0 

-W  r  dBze (r, t) 

+  a25  0 

0 
0 

Now,  by  using  Equation  (3.2)  in  Equation  (3.3),  the 
following  five  differential  equations  result. 


<  d/dt  )  *N(r,t) 


<3. 4) 


• •--»  v.— tr* y*  w  -  > ' 


=  ZN°_  $Vr(r,t)  +  r  Dr  (tfVr(r,t>) 


r  u 

<  d/dt  )  <fVr(r,t)  =  a23  SVq( r,t)  +  934  r  $Vz(r,t>  <3. 5) 

+  a25  *Er(r,t)  -  W  r  $Bze(r,t> 

<  d/dt  >  $Ve( r,t)  =  -a23  $Vr( r,t)  (3.6) 

<  d/dt  )  *Vz(r,t)  =  -a24  r  tfVr(r,t)  (3.7> 

(  d/dt  )  $Er(r,t)  =  a52  $Vr(r,t)  C3.8) 


3-3  The  Numerical  Approx i mat i on . 

Equations  <3. 4)  through  (3.8)  are  the  linear,  time 
invariant,  differential  equations  that  govern  the  dynamic 
behavior  of  the  particle  beam.  However,  only  Equation  <3. 4) 
contains  a  partial  derivitive  with  respect  to  r.  This 
partial  derivitive  is  replaced  with  its  finite  difference 
approximation  (backward  differences)  prior  to  investigating 
the  solutions  to  Equations  (3.4)  through  (3.8).  This  finite 
difference  approximation  is  as  follows  : 


Dr(*Vr(r,t))  *  (*Vr(r,t)  -  $Vr(r-h,t>)  /  h 


(3.9) 


Where  h  is  the  spacing  between  points.  Therefore,  Equation 
(3.4)  becomes: 


<d/dt)  *N(r,t)=  , 

<-N°  /r)|VVr(r,t>  +<r/h)  <  *Vr(r,t)  -  $Vr(r-h,t)  )J  <3.1i 

3-4  Solution  Technique. 


Out  of  the  five  governing  differential  equations 
discussed  earlier,  only  Equation  (3.10)  requires  a  finite 
difference  approximation.  Now,  since  none  of  the  other  four 
equations  are  coupled  to  Equation  (3.10),  their  solutions 
can  be  found  by  applying  standard  differential  equation 
techniques.  The  solutions  to  these  four  equations  are  as 
follows  : 


*Vr(r,t>  = 


*V0(r,t)  = 


$V,(r,t)  * 


6Er(r,t)  = 


a25  a52  r  sln[  Qt  1 
100  (2 


a23  a25  a52  r  <  cost  Qt  1  -  1  > 

100  a 


a24  a25  a52  r2  <  cost  Qt  1  -  1  > 
100  Q 


-  a522  a25  r  (  cost  Qt  1  -  1  ) 
100  Q 


(3.11) 


(3. 12) 


(3.13) 


(3.14) 
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where 


a  =  <  a232  +  ^242  <r2)_  352  *25  >1/2 

Note  that  this  yields  the  same  values  for  ft  as  does 
Equat ion  (2.8). 

It  should  be  noted  that  in  the  solution  of  these 
four  equations,  the  initial  values  are  zero.  This  is 
consistent  with  the  derivation  of  the  perturbation  equations 
and  agrees  with  the  closed  form  solution  found  in  Chapter  2 
for  t=0. 

The  solution  of  Equation  (3.10),  however,  is  slightly 
more  difficult.  Since  Equation  (3.10)  contains  finite 
differences  and  is  coupled  to  the  other  four  equations,  the 
equation  must  be  analyzed  at  specific  points  along  the 
radius  of  the  beam. 

Now,  after  substituting  in  Equation  (3.11),  Equation 
(3. 10)  becomes 

-N0  f  a25  a52  r  sint  ftt  1 
(  d/dt  )  *N(r,t)  =  -  - 

r  100  ft 

1  — 

r  a25  a52  r  sini  ftt  1  (3.15) 

+  . .  . . . 

100  h  ft 

a25  a52  <r-h>  slnl  ftl  t  ] 

100  h  ftl 

where 
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..  . - .  - ,  - . - k\;. ,\ 


Q1  =  C  a232  +  a242  <r  -  h)2  -852  a25  ) 


Now  the  solution  of  Equation  (3.15)  can  be  found 
through  simple  integration,  where  the  resulting  constant  of 
integration  is  determined  by  using  the  fact  that  dN(r,t>  is 
zero  at  t=0.  Therefore,  the  solution  of  Equation  (3.15)  is: 


$N(r,t)  = 


-N°  fa25  a52  <  1-cos C  Ot  1  ) 


(3. 16) 


100  0 


r  < 1-cosC  Qt  1  ) 


h  01 


<r-h)  <  1-cost  01  t  1  > 


h  01 


We  now  have  the  solutions  to  all  five  governing  differential 
equations. 


3-5  Data  Presentation. 


In  order  to  compare  the  results  from  Chapter  3  with 
those  found  in  Chapter  2,  the  radius  of  the  partlcal  beam  is 
divided  into  twenty  evenly  spaced  intervals.  Since  the  beam 
radius  was  given  as  ten  millimeters,  the  spacing  (  h  )  is 
one  half  of  a  millimeter.  For  each  value  of  time  chosen 
there  will  be  nineteen  results  for  the  five  equations 
derived  in  Section  3-4.  Each  of  these  results  represents  a 


* 


p-VS 

4 


V*  '»*.  **. 
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solution  at  a  specific  radial  distance  from  the  center  of 
the  beam.  The  same  format  will  be  used  for  the  data 
presented  from  Chapter  2. 

The  results  from  this  chapter  are  presented  in  Table 
thru  20  and  are  shown  graphically  in  Figures  14  thru  17. 


CHAPTER  4 


4- 1  Introduct Ion 

The  material  presented  in  this  chapter  is  an  expansion 
of  that  presented  in  Chapter  3  for  the  system  (Equation  2.1) 
which  results  when  the  Electrostatic  Approximation  Model  is 
not  used  to  reduce  the  complexity  of  the  system.  The 
resulting  ninth-order  system  is  solved  and  compared  with  the 
results  of  Chapter  3  in  order  to  evaluate  the  effect  of  the 
Electrostatic  Approximation. 


4-2  Setting  up  the  Problem 

The  ninth-order  system  of  governing  differential  equations 
are 


<d/dt>  X(t>  =  A  X(t)  +  g< t)  <4. 

where  A  and  g(t>  are  defined  in  Chapter  2,  Equations 
(2.1b)  and  (2.1c>. 

As  In  Chapter  3,  the  only  external  field  applied  is 
<fBz®(r,t),  and  its  maximum  value  is  assumed  to  be  one 


percent  of  the  equilibrium  value  B0,  which  was  previously 
determined  in  Chapter  2.  The  rest  of  the  external  fields 
are  assumed  equal  to  zero.  Therefore,  from  Equation  <2. lb) 


g<  t )  = 


0 

-  W  r  <0.01)  B0 
0 
0 
0 
0 
0 
0 
0 


Multiplying  out  and  combining  the  terms  described  in 
Equation  (4.1)  results  in  the  following  nine  partial 
differential  equations: 


<d/dt) 

*N<r,t) 

=  -N°  Dr(tfVr<r,t> ) 

(4.2) 

.  < 

■  %  w  ’ 

<d/dt) 

t' Vr  <  r ,  t ) 

+  a25 

-  a25 

=  a23  $V9(r,t)  +  a24  r  6  Vz(r,t> 

*Er<r,t)  -  a25  Vz°  $B9<r,t) 

W  r  6Bz(r, t)  +  W  r  (.01)  B0 

(4.3) 

Sv-v 

‘  *  *’» 

'*/ 

v  «o  *<, 

<d/dt> 

*V9< r,t) 

=  -a23  #Vr<r,t)  +  a25  £E9(r,t) 

(4.4) 

y*\V 

aV/ 

If- .  JL 

<d/dt) 

*Vz<r,t> 

=  -a24  r  JVr<r,t>  +  a25  £Ez(r,t) 

(4.5) 

where  h  Is  the  spacing  between  points.  After  Implementing 
this  finite  difference  scheme,  the  five  equations  discussed 
above  result  in  the  following: 


<d/dt) 


<d/dt ) 


(d/dt) 


(d/dt) 


<d/dt) 


$N(r,t)  =  ^N0 

r 


+  r 


$Vr(r,  t> 

tfVr(r,t)  -  $Vr( r,t> 


(4.11) 


JEe<r,t)  =  -a52  W  r  $N(r,t)  +  a52  *Ve(r,t)  (4.12) 


N° 


-  C2  tfB,(r,t)  -  tfB,(r-h,t) 


$Ez(r,t)  =  a52  Vz°  tfN < r , t )  +  a52  $Vz(r,t) 
N° 


(4. 13) 


+  C2  tfB0(r,t)  +  r 
r 


«SB0(r,t)  -  dB0(r-h,  t) 


tfB0(r,t)  =  6Ez(r,t)  -  <9Ez(r  -  h,t) 

h 


(4.14) 


*Bz(r,t)  = 


-J,  4E0(r,t)  +  rr#E0(r,t)  -  dE©(r  -  h,t>~| 


(4.15) 


4-4  Final  Analysis  and  Solution  Techni 


We  now  have  at  hand,  and  In  the  correct  form,  the  nine 


simultaneous  differential  equations  that  we  will  solve. 


Since  these  equations  are  much  more  complicated  than  those 
solved  in  Chapter  3,  a  numerical  integration  routine  was 
used.  The  results  from  this  integration  routine  are 
tabulated  numerically  in  Tables  21  and  22.  Figure  20  shows 
how  the  tfN(r,t)  diverges  with  time.  Since  all  nine 
solutions  diverge  in  a  similar  manner  no  other  figures  are 
presented. 


CHAPTER  5 


RESULTS  AND  CONCLUSIONS 


5 . 1  Introduct  ton. 


This  chapter  contains  a  discussion  of  the  results 
obtained  from  the  previous  four  chapters.  Based  upon  these 
results,  conclusions  are  drawn  and  recommendations  are 
presented . 


5-2  Closed  Form  Solution. 


5-2.1  Step  Input. 

The  results  from  the  analysis  presented  in  Section  2-4 
are  tabulated  in  Tables  1  thru  4  and  are  shown  graphically 
in  Figures  2  thru  5.  It  is  observed  from  Figure  2  that  a 
stable  solution  for  *N(r,t>  is  realized  with  no  significant 
damping  apparent.  Identical  results  are  seen  from  the 
equations  which  provide  the  solutions  to  $Vz(r,t)  and 
£Er(r,t)  (Figures  4  and  5).  This  is  consistent  with  the 
form  of  the  solutions  found  in  Equation  (2.11).  As  can  be 
seen  in  Equation  (2.11)  the  solutions  £N(r,t),  $Vz(r,t),  and 
$Er(r,t)  all  contain  the  same  trigonometric  function,  (1-cos 


tfltl).  The  only  difference  in  the  solutions  is  the 
amplitude  of  the  trigonometric  term;  therefore,  similar 
results  are  expected.  Figure  3  shows  that  the  solution  of 
tfVr(r,t)  is  also  stable  with  no  damping  apparent. 


An  important  result  of  this  section  is  that  the  closed 
form  solutions,  as  determined  in  Reference  tli,  can  be  used 
to  evaluate  the  solutions  to  the  system  in  question,  and 
that  these  results  make  physical  sense. 

5-2.2  Input  Rapidly  Increasing  to  Constant. 

The  results  from  the  analysis  presented  in  Section 
2-5.1  are  tabulated  in  Tables  5  thru  8,  and  are  shown 
graphically  in  Figures  6  thru  9.  It  should  be  noted  that 
this  type  of  input  (Figure  18)  is  similar  to  the  step  input 
of  Section  2-4.  These  two  types  of  inputs  are  unequal  only 
over  a  relatively  short  interval  of  time,  approximately  from 
t=0 . 0  to  t=4.0xl0-9  seconds,  for  the  number  of  significant 
digits  used  in  this  effort.  Therefore,  it  makes  physical 
sense  that  the  results  of  Section  2-5.1  are  similar  to  those 
of  Section  2-4.  The  lag  in  the  solution  that  is  apparent 
when  comparing  the  data  from  Section  2-4  to  the  data  of 
Section  2-5.1  is  conjectured  to  be  due  to  a  slight  lag  in 
the  values  of  the  input  due  to  the  exponential  term.  Also, 
Figures  6  thru  9  show  that  the  solutions  from  this  section 
are  stable  and  exhibit  no  apparent  damping. 


5-2.3  Input  Rapidly  Increasing  to  Constant  Followed  by  a 


Rapidly  Decrease  to  Zero. 

The  results  from  the  analysis  presented  In  Section 
2-5.2  are  tabulated  in  Tables  9  thru  16  and  are  shown 
graphically  In  Figures  10  thru  13.  These  figures  clearly 
show  two  distinct  functions  for  the  solutions,  X<t) .  The 
first  function  in  all  four  figures  spans  the  time  Interval 
between  t=0.0  and  t=4.5xl0-7  seconds.  Over  this  Interval, 
the  solutions,  X(t),  are  almost  Identical  to  those  values 
found  in  Section  2-5.1.  This  can  be  seen  by  comparing  the 
values  in  Tables  9,11,13,  and  15  with  their  related  values 
in  Tables  5  thru  8.  Actually,  the  two  solutions  for  X(t) 
are  exact  in  comparison  between  t=0.0  and  t  =  4. 45xl0“7 
seconds,  with  the  only  differences  encountered  between  t  > 
4.45xl0~7  to  t=4.5xl07  seconds.  This  can  be  explained  by 
referring  to  Equation  <2.17),  which  shows  how  the  external 
axial  magnetic  input  was  defined.  Since,  the  input  is 
exactly  the  same  as  that  used  for  Section  2-5.1  over  the 
interval  0  <  t  <  tj,  the  results  are  expected  to  be 
identical.  Over  the  interval  tj  <  t  <  t2  the  input  used  is 
the  cubic  function  defined  in  Equation  (2.17b).  The 
solutions  are  found  by  summing  the  results  of  Equation  (2.6) 
integrated  from  t=0 . 0  to  t=tj,  where  the  input  is  defined  by 
Equation  (2.17a),  with  the  results  of  Equation  (2.6) 
integrated  from  t=tj  to  t<t2>  where  the  input  is  defined  by 


integrated  from  t=tj  to  t<t2»  where  the  input  is  defined  by 
Equation  <2. 17b).  Therefore,  the  results  over  this  interval 
are  expected  to  be  slightly  different.  The  second  function 
in  all  four  figures  spans  the  interval  of  t  >  4.5xl0~7  to 
t=«  seconds.  The  input  used  over  this  interval  is  equal  to 
zero;  therefore  the  solutions  are  the  summation  of  the 
previous  works  integrated  over  their  associated  intervals. 
Both  these  functions  encountered  on  Figures  10  thru  13  are 
stable  and  exhibit  no  apparent  damping. 

5-3  Numerical  Solution  Technique. 


The  results  from  the  analysis  presented  in  Chapter  3 
are  tabulated  in  Tables  17  thru  20  and  are  shown  graphically 
in  Figures  14  thru  17.  From  these  figures  it  is  concluded 
that  the  solutions,  X<t),  are  stable  and  show  no  apparent 
damping.  It  can  also  be  concluded  by  comparing  the  data  in 
Tables  17  thru  20  with  their  related  data  from  Section  2-4, 
that  the  closed  form  solution  techniques  derived  in 
Reference  ill  produces  excellent  results.  This  conclusion 
is  verified  by  the  results  of  Chapter  3.  It  is  conjectured 
that  the  small  differences  found  between  the  closed  form 
solution  technique  of  Section  2-4  and  the  numerical  solution 
technique  of  Chapter  3  are  due  to  errors  associated  with  the 
Chapter  3  method.  Under  this  assumption,  if  the  number  of 
points  across  the  radius  of  the  beam  were  increased  in  order 


to  increase  the  accuracy  of  this  solution  technique,  the 
results  of  this  would  more  closely  approximate  those  of 
Section  2-4. 


5-4  Numerical  Solution  Technique  Without  Electrostatic 
Approx i mat  1  on 

The  results  of  the  Chapter  4  analysis  is  presented  in 
Tables  21  and  22.  The  solution  for  SN(r,t>  is  shown 
graphically  in  Figure  20.  It  is  observed  from  Figure  20  and 
from  the  data  presented  in  Tables  21  and  22  that  the 
solutions,  X(t),  are  not  stable  for  all  time  t.  However, 
the  data  in  Tables  21  and  22  do  approximately  agree  with 
that  of  Section  2-4  over  a  Interval  of  time,  t=0.0  to 
t<3.0xl0-11  second.  Over  this  interval  it  can  be  concluded, 
due  to  the  agreement  in  associated  values,  that  the 
Electrostatic  Approximation  is  an  adequate  one.  For  times 
other  than  that  described  above,  the  solution,  X(t>, 
diverges  to  t  •.  It  is  conjectured  that  this  instability  is 
being  caused  by  any  of  the  following  anomalies.  First,  a 
numerical  instability  could  have  been  encountered  by  the 
integration  routine  used  to  solve  the  differential 
equations.  Secondly,  the  process  of  linearizing  the 
differential  equations  might  have  eliminated  a 
damping/restoring  term  present  in  the  nonlinear  equations. 
This  term  would  had  to  have  been  in  either  £EQ(r,t), 


£Ez<r,t),  £Br(r,t),  £B0(r,t),  or  $Bz<r,t),  since  the  Section 
2-4  material  was  not  effected  in  a  similar  manner.  However, 
it  is  more  likely  that  in  choosing  B0  such  that  we  have  the 
the  minimun  external  magnetic  field  required  to  maintain 
rigid  rotor  equilibrium,  we  produce  a  system  which  borders 
on  Instability.  Therefore,  an  increase  in  B0  would  tend  to 
stabilize  this  system. 


A 


5-5  Recommendat ions 

The  following  recommendations  are  made  as  a  result  of 
this  effort. 

1.  Investigate  the  divergence  of  the  Chapter  4  study. 

1.  reexamine  the  derivation  of  the  dynamic  model 

found  in  Reference  Cil  paying  close  attention  to 
the  last  four  terms. 

ii.  continue  work  of  chapter  4  by  increasing  B0  and 
recalculating  the  results.  Should  result  in 
stable  solutions. 

2.  Investigate  a  closed  form  solution  to  the  dynamic  system 
found  in  Reference  (11  proir  to  the  implementation  of 
the  electrostatic  approximation. 
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TABLE  1 

PERTURBED  NUMBER  DENSITY,  SECTION  2-4 


T 1 

T2 

T3 

T  4 

T5 

R 

7.668580*12 

1 .782200*13 

1 .344340*13 

l .871190*12 

2.500090*12 

5 . 00000E-04 

7.668580*12 

1 .782200*13 

1 . 344340*13 

1.871190+12 

2 . 50009D+ 1 2 

1 .00000E-03 

7.668580*12 

1.782200*13 

1 . 34434D* l 3 

1 .871 190+12 

2.500090*12 

1 .50000E-03 

7.668580*12 

1.782200*13 

1  . 344340*13 

1.871190*12 

2.500090*12 

2 . OOOOOE-03 

7.668580*12 

1 .782200*13 

1  . 344340*13 

1 .871 19D+ 12 

2.500090*12 

2.50000E-03 

7.668580*12 

1.782200*13 

1  . 34434D+13 

1.871 190*12 

2 . 50009D+ 1 2 

3. OOOOOE-03 

7.668580*12 

l .782200*13 

1 . 344340*13 

1 .871190*12 

2.500090+12 

3 . 50000E-03 

7.668580*12 

1.782200*13 

1 . 344340*13 

1.871  1 90*  1  2 

2 . 500090+ 1 2 

4 . OOOOOE-03 

7.668580*12 

1 .782200*13 

1 .344340*13 

1 .871 19D*12 

2 .500090+ 1 2 

4 . 50000E-03 

7.668580*12 

1.782200*13 

1.344340*13 

1.871190+12 

2 . 500090+ 1 2 

5. OOOOOE-03 

7 . 66858D* 1 2 

1 .782200*13 

1 . 34434D* 1 3 

1.8711 90* 1 2 

2.500090+12 

5 . 50000E-03 

7.668580*12 

1 .78220D*13 

1.344340*13 

1 .871 19D+12 

2.500090+12 

6 . 00000E-03 

7.668580*12 

1 .782200*13 

1.344340*13 

1 .871190*12 

2 . 50009D+ 1 2 

6.S0000E-03 

7.668580*12 

1.782200*13 

1.344340*13 

1.87U9D+12 

2 . 50009D* 1 2 

7. OOOOOE-03 

7.668580*12 

1.782200*13 

1.344340*13 

1.871190*12 

2.500090*12 

7.50000E-03 

7.668580*12 

1.782200*13 

1.344340*13 

1.871190*12 

2 . 50009D+ 1 2 

8. OOOOOE-03 

7.668580*12 

1.782200*13 

1 .344340*13 

1.8711 90* 1 2 

2 . 50009D* 1 2 

8.50000E-03 

7.668580*12 

1 ,78220D*13 

1.344340*13 

1.8711 90* 1 2 

2 . 50009D+ 1 2 

9. 00000E-03 

7 . 66858D+  t  2 

1.782200*13 

l . 344340*13 

1.871190*12 

2.500090*12 

9 . 50000E-03 

TABLE  2 

PERTURBED  RADIAL  VELOCITY,  SECTION  2-4 


Ti  T2  T3  7  4  T5  R 


-  1 . 38946D+02 
-2.778930+02 
-4. 168390+02 
-3 . 557860+02 
-6.947320+02 
-8.336790+02 
-9.726250+02 
-1.11 1570+03 
-1.250520+03 
-1 . 389460+03 
-1 .52841D+03 
-I .667360+03 
-1 .80630D+03 
-1 .94323D+03 
-2.084200+03 
-2 . 2231 4D+03 
-2. 36209D+03 
-2 . 30 104D+03 
-2 . 63998D+03 


-4.502240*01 
-9.004480+01 
-1.350670+02 
-1.800900+02 
-2.251120+02 
-2.701340+02 
-3. 15 157D+02 
-  3 . 60 1 79D+02 
-4.03202D+02 
-4.302240+02 
-4.952470+02 
-3.402690+02 
-3.852910+02 
-6.303140+02 
-6.73336D+02 
-7.203590+02 
-7.653810+02 
-8. 104030+02 
-8.354260+02 


1 . 243580+02 
2.487160+02 
3.73074D+02 
4.974320+02 
6.217900+02 
7.461480+02 
8.703060+02 
9.948640+02 
1 . 119220+03 
1 . 243580+03 
1.367940+03 
1.492300+03 
1.616650+03 
1.741010+03 
1.863370+03 
1.989730+03 
2.  1  1409D+03 
2.238440+03 
2.362800+03 


8 . 53 l 78D+0 1 
1 . 70636D  +  02 
2.539530+02 
3.412710+02 
4.263890+02 
3.  119070+02 
5 . 97224D+02 
6 . 823420+02 
7 . 678600+02 
8.331780+02 
9 . 38493D+02 
1 .02381D+03 
1 . 109130+03 
1 . 194450+03 
1 .27977D+03 
1.363080+03 
1 .43040D+03 
1.333720+03 
1.621040+03 


-9.671 27D+0 1 
-1 . 93425D+02 
-2.901380+02 
-7  8683 10+02 
-4.833640+02 
-5 . 80276D+02 
-6.769890+02 
-7 . 73702D+02 
-8.704130+02 
-9.671270+02 
- 1 . 06384D+03 
-1 . 160330+03 
-I. 237270+03 
- 1 . 353980+03 
-1.450690+03 
-1.347400+03 
-  1 . 644 1 2D+03 
-1.740830+03 
-1 .83734D+03 


3 . 00000E-04 
1 . 00000E-03 
1 . 50000E-03 
2 . 00000E-03 
2 . 50000E-03 
3 . 00000E-03 
3 . 50000E-03 
4 . OOOOOE-03 
4 . 50000E-03 
3. 00000E-03 
3 . 30000E-03 
6. OOOOOE-03 
6 . 30000E-03 
7. 00000E-03 
7 , 30000E-03 
8. 00000E-03 
8 . 50000E-03 
9. OOOOOE-03 
9 . 30000E-03 


TABLE  3 

PERTURBED  AXIAL  VELOCITY,  SECTION  2-4 


4 . 0  t  3350-04 
I .60534D-03 
3.612020-03 
6.421370-03 
1.003340-02 
1 .444810-02 
1 . 966340-02 
2 . 56835D-02 
3.250820-02 
4.01 335D-02 
4.836160-02 
5 . 77923D-02 
6.782370-02 
7 . 866 1 7D-02 
9.030050-02 
1.027420-01 
1 . 139860-01 
1 . 300330-01 
I . 448820-01 


-9. 327140-04 
-3.730860-03 
-8.394430-03 
-1.492340-02 
-2 . 331790-02 
-3.337770-02 
-4.370300-02 
-5. 969370-02 
-7.554990-02 
-9 . 327  1  40-02 
-1 . 1 28380-0 1 
-1.343110-01 
-1 .37629D-01 
-1 .828120-01 
-2.098610-0' 
-2.387750-01 
-2.693540-01 
-3.021990-01 
-3.367100-01 


-7.035600-04 
-2.814240-03 
-6 . 332040-03 
-1 . 123700-02 
-1 .738900-02 
-2 . 53282D-02 
-3.447440-02 
-4.30278D-02 
-3.698840-02 
-7.033600-02 
-8.313080-02 
-1.0131 3D-01 
-1 . 18902D-01 
-1 . 378980-01 
.583010-01 
-  1.8011 10-01 
-2.033290-01 
-2.279530-01 
-2.339830-01 


-9 . 79289D-05 
-3.9171 6D-04 
-8.81 3600-04 
-  1 . 366860-03 
-2 . 44822D-03 
-3.32344D-03 
-4.798520-03 
-6.267430-03 
-7.932240-03 
-9 . 79289D-03 
-1 . 184940-02 
-1 . 41018D-02 
-1.633000-02 
-1.919410-02 
-2.203400-02 
-2.306980-02 
-2.830150-02 
-3. 17290D-02 
-3.533230-02 


-1 . 308420-04 
-5 . 23369D-04 
-1 . 177380-03 
-2.093480-03 
-3.271 06D-03 
-4 . 7 1032D-03 
-6.41 127D-03 
-8 . 3739  ID-03 
-1.059820-02 
-1 . 308420-02 
-1 .S8319D-02 
-1.884130-02 
-2.211240-02 
-2.564510-02 
-2.943930-02 
-3.349360-02 
-3.781340-02 
-4.239290-02 
-4.723410-02 


5.00000E-04 
1 . 00000E-03 
1 . 30000E-03 
2 . OOOOOE-03 
2 . 30000E-03 
3. 00000E-03 
3 . 50000E-03 
4. OOOOOE-03 
4 . 30000E-03 
3. OOOOOE-03 
3 . 30000E-03 
6. 00000E-03 
6 . 30000E-03 
7. 00000E-03 
7 . 30000E-03 
8. OOOOOE-03 
8 . 50000E-03 
9. 00000E-03 
9 . 50000E-03 


TABLE  4 

PERTURBED  RADIAL  ELECTRIC  FIELD,  SECTION  2-4 


T 1 

T2 

T3 

T  4 

T5 

R 

3. 469070*01 

8.062220*01 

6.081450*01 

8 . 4648 1D+00 

1 . 130980+01 

5 . 00000E-04 

6 .938150*01 

1.612440+02 

1.216290*02 

1 .692960*01 

2.261960+01 

1 . 00000E-03 

1.040720*02 

2.418670*02 

1.824430+02 

2 . 3  39440  +  01 

3.392930+01 

1 . 50000E-03 

1 . 387630*02 

3.224890+02 

2.43258D+02 

3 ■ 38392D+0 1 

4.523910+01 

2 . 00000E-03 

1 .734540*02 

4.0311 10*02 

3.040720*02 

4 . 232400+0 1 

3 . 65489D+0I 

2.30000E-03 

2.081440*02 

4.837330*02 

3.648870+02 

5.078880+01 

6 . 783870+0! 

3 . 00000E-03 

2.428350*02 

5.643550+02 

4. 257010*02 

5 . 925360+0 1 

7.916830+01 

3 . 50000E-03 

2.775260*02 

6.449780*02 

4.863160+02 

6.771850+01 

9.047830+01 

4.00000E-03 

3. 122170*02 

7.256000*02 

5.473300*02 

7.61833D+OI 

1 . 01788D+02 

4.30000E-03 

3.469070+02 

8.062220+02 

6.081430+02 

8 . 464810+01 

1 . 1 3098D+02 

3.00000E-03 

3.815980*02 

8.868440+02 

6.689390+02 

9.31 129D+01 

1 . 24408D+02 

5 . 50000E-03 

4. 162890*02 

9 . 67467D+02 

7.297740+02 

1 .013780+02 

1 . 33717D+02 

6 . 00000E-03 

4. 50980D*02 

1.048090+03 

7.903880+02 

1 . 10042D+02 

1.470270+02 

6 . 30000E-03 

4 . 85670D+02 

1 . 12871D+03 

8.314030+02 

1 . 183070+02 

1.383370+02 

7 . OOOOOE-03 

5.203610*02 

1 -20933D+03 

9.122170+02 

1.269720+02 

1.696470+02 

7 . 50000E-03 

5.530520*02 

l .289960+03 

9.730320+02 

1.334370+02 

1 . 80957D+02 

8. OOOOOE-03 

5.897420*02 

1 . 37058D+03 

1 .033850+03 

1 . 439020+02 

1 .92266D+02 

8 . 50000E-03 

6.244330+02 

1.451200+03 

1 .094660+03 

1.523670+02 

2.035760+02 

9. 00000E-03 

6.591240+02 

1 .53182D+03 

1 . 155480+03 

l .608310+02 

2.148860+02 

9 . 50000E-03 

TABLE  5 

PERTURBED  NUMBER  DENSITY,  SECTION  2-5 


7.413240*12 

1.774640*13 

l . 367410*13 

2.021430+12 

2.31 793D+ l 2 

5.00000E-04 

7.41 324D  + I 2 

1 .774640*13 

1 . 367410* 1 3 

2.021430+12 

2 . 3 1 793D+  I  2 

1 .OOOOOE-03 

7.413240*12 

1.774640*13 

1 .367410+13 

2.021 43D+ 1 2 

2 . 3 1 79  3D+  l  2 

1 . 50000E-03 

7.413240*12 

1.774640*13 

1.367410+13 

2.021430+12 

2.31 793D  + l 2 

2. OOOOOE-03 

7.413240*12 

1.774640*13 

1.367410+13 

2.021 43D+ 1 2 

2.317930+12 

2.50000E-03 

7.41 324D+ 1 2 

1.774640+13 

1.367410+13 

2.021 43D+ 1 2 

2.317930+12 

3. OOOOOE-03 

7.41 324D* 1 2 

1.774640*13 

1.367410+13 

2.021430+12 

2.3I793D+12 

3.50000E-03 

7.413240*12 

1 .77464D+13 

1.367410*13 

2.021 43D+ 1 2 

2.317930+12 

4. OOOOOE-03 

7.413240*12 

1.774640+13 

1.367410+13 

2.021430+12 

2. 31793D+12 

4.30000E-03 

7.413240*12 

1.774640+13 

1.367410+13 

2.021 43D+ 1 2 

2.317930+12 

3. OOOOOE-03 

7.413240*12 

1 .77464D+13 

1.367410+13 

2.021430+12 

2.31 793D+ l 2 

5.S0000E-03 

7.413240*12 

1.774640+13 

1 . 3674  ID*  1  3 

2.021430+12 

2.317930+12 

6. 00000E-03 

7.413240*12 

1 .774640+13 

1.367410+13 

2.021430+12 

2. 31793D+12 

6.50000E-03 

7.413240+12 

1 .77464D+13 

1 .36741D+13 

2.021430+12 

2 . 31793D+ 1 2 

7. OOOOOE-03 

7.413240+12 

1.774640+13 

1.367410+13 

2.021430+12 

2.31793D+12 

7 . 50000E-03 

7.413240+12 

1.774640+13 

1 . 3674 1D+ 1 3 

2.021430+12 

2.317930+12 

8. 00000E-03 

7.413240+12 

1 .77464DH3 

1.367410+13 

2. 02 t 430+ 12 

2 . 31793D+ l 2 

6 . 50000E-03 

7.413240+12 

l .77464D+13 

1 .36741D+13 

2.021430+12 

2.317930+12 

9. OOOOOE-03 

7.413240+12 

1.774640*13 

1.367410+13 

2.021430+12 

2.317930+12 

9 . 50000E-03 

TABLE  6 

PERTURBED  RADIAL  VELOCITY ,  SECTION  2-B.1 


Tl 

T2 

T3 

1 . 38 193D*02 

-4 . 87404D+01 

1 . 224000+02 

2 . 763860*02 

-9.748090+01 

2 .448000+02 

4. 1 45790+02 

-l  .  4622 ID  +  02 

3.672000+02 

5  -52773D+02 

-I  .949620*02 

4.896000+02 

6. 909660+02 

-2.437020+02 

6. 120000+02 

8  .  291590+02 

-2.924430*02 

7.344000+02 

9 . 67352D+02 

-3.41 1830+02 

8 . 968000+02 

1  .  105550+03 

-3 . 89924D+02 

9.792000+02 

1  . 24374D+03 

-4 . 38664D+02 

1 . 101600+03 

t . 381930+03 

-4.874040*02 

1.224000+03 

1  -52013D+03 

-5.361 45D*02 

1.346400+03 

l  . 658320+03 

-5.848850*02 

1.468800+03 

1 .79651D+03 

-6 . 336  26D*02 

1.591200+03 

t . 934700+03 

-6.823660*02 

1 .713600+03 

2.072900*03 

-7.311 07D*02 

1 .83600D+03 

2.211 09D*03 

-7.798470*02 

1.938400+03 

2 . 34928D+03 

-8.285880*02 

2 . 08080D+03 

2 . 48748D+03 

-8.773280*02 

2.203200+03 

T  4 

T5 

R 

8.8401 3D+0 1 

-9. 373550+01 

5 . 00000E-04 

1 .768030  +  02 

-1.8731 10+02 

1 . OOOOOE-03 

2.652040*02 

-2.812670+02 

1 .50000E-03 

3.536030+02 

- ? . 750220+02 

2. 00000E-03 

4.420070+02 

-4 . 68778D+02 

2 . 50000E-03 

5 . 304080  +  02 

-3.623330+02 

3. 00000E-03 

6.188090+02 

-6 . 56289D+02 

3. 50000E-03 

7.072110+02 

-7.500440+02 

4. OOOOOE-03 

7.936120+02 

-8.438000+02 

4 . 30000E-03 

8.840130+02 

-9. 37533D+02 

5. OOOOOE-03 

9.724150+02 

-1 .031310+03 

5 . 30000E-03 

1 . 06082D+03 

-1.123070+03 

6 . 00000E-03 

1 . 1 49220+03 

-1.218820+03 

6 . 50000E-03 

1.237620*03 

-1.312580+03 

7. 00000E-03 

1.326020+03 

-1 . 40633D+03 

7 . 50000E-03 

1.414420+03 

-1 .30009D+03 

8. 00000E-03 

1 .302820+03 

-1.593840+03 

8 . 30000E-03 

1.591220+03 

-1.687600+03 

9. OOOOOE-03 

2.623670*03-9.260680+02  2.32560D+03  1.679630*03-1.781330*03  9.30000E-03 


T'ASLE  7 

PERTURBED  .AXIAL  VELOCITY,  SECTION  2-5.1 


T 1  T2 

-3. 879720-04  -9.287560-04 
-1.551890-03  -3.715020-03 
-3.491750-03  -8.358800-03 
-6.207550-03  -1.486010-02 
-9.699290-03  -2.321890-02 
-1.396700-02  -3.343520-02 

-  1 . 90 1 06D-02  -4.550910-02 
-2.483020-02  -5.944040-02 
-3.142570-02  -7.522930-02 
-3.879720-02  -9.287560-02 
-4.694460-02  -1.123800-01 
-5.586790-02  -1.337410-01 
-6.556720-02  -1.569600-01 
-7.604250-02  -  1 . 82036D-0 1 
-8.729370-02  -2.08970D-01 
-9.932080-02  - 2 . 3776 2D-0 1 
-1.121240-01  -2.684110-01 
-1.257030-01  -3.009170-01 

-  1 . 40058D-0 1  -3 . 3528 10-0  t 


T3  T4 

-7.156330-04  -1.057920-04 
-2.862530-03  -4.231670-04 
-6.440700-03  -9 . 52 1 26D-04 
-1. 145010-02  -  1 . 69267D-03 
-1.789080-02  -2 . 64479D-03 
-2.576280-02  -3.808500-03 
-3.506600-02  -5.18380D-03 
-4.580050-02  -6.770670-03 
-5.796630-02  -8.569130-03 
-7 . 1 56  340-02  -  1 . 05792D-02 
-8.659170-02  -1.280080-02 
-1.030510-01  - 1 . 52340D-02 
-1.209420-01  -1.787880-02 

-  1 . 40264D-01  -2 . 07352D-02 
-1.810180-01  -2.380310-02 

-  I . 83202D-0 1  -2.708270-02 
-2.068180-01  -3.057380-02 
-2.318650-01  -3.427650-02 
-2.583440-01  -3.81908D-02 


T5  R 

-1.213090-04  5 . 0O000E-04 
-4.852360-04  1.00000E-03 
-1.091780-03  1 . 50000E-03 
-  1 . 94094D-03  2.00000E-03 
-3.032730-03  2.50000E-03 
-4.367120-03  3 . 00000E-03 
-5.9441 4D-03  3.50000E-03 
-7 . 76378D-03  4.00000E-03 
-9 . 826030-03  4.50000E-03 
-1.213090-02  5 . OOOOOE-03 
-1.467840-02  5 . 50000E-03 
-1.746850-02  6. OOOOOE-03 
-2 . 0501 2D-02  6 . 50000E-03 
-2.377660-02  7. 00000E-03 
-2 . 72945D-02  7.50000E-03 
-3. 1055  ID-02  8. 00000E-03 
-3.505830-02  8.50000E-03 
-3 . 9304  ID-02  9. 00000E-03 
-4 . 37926D-02  9.S0000E-03 
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PERTURBED  RADIAL  ELECTRIC  FIELD,  SECTION  2-5.1 


T! 

T2 

T3 

T  4 

T5 

R 

3 . 35  3560*0 1 

8.028010*01 

6 . 1 858 1D+01 

9. 144450*00 

1 . 04857D+0 I 

5 . 00000E-04 

6.707120+01 

1 .603600*02 

I .237160*02 

l .828890*01 

2.0971 5D+0 1 

l . 00000E-03 

1 . 006070+02 

2.408400*02 

1 . 85574D+02 

2.743340*01 

3. 145720*01 

1 . 50000E-03 

l  . 341420*02 

3.21 1200*02 

2.474320*02 

3.657780*01 

4. 194300+01 

2. 00000E-03 

I  .  676780*02 

4.014000*02 

3.092910*02 

4 . 57223D+0 1 

3 . 24287D+01 

2.50000E-03 

2.0121 40+02 

4.816800*02 

3.71 1490*02 

3 . 48667D+0 1 

6 . 291440+01 

3.00000E-03 

2 . 347490*02 

3.619600+02 

4 . 330070+02 

6 . 401120+01 

7.340020+01 

3 . 50000E-03 

2 . 68285D+02 

6.422410+02 

4 . 94863D+02 

7 . 3 1336D*0 l 

8.388390+01 

4.00000E-03 

3.018200*02 

7.225210+02 

3.367230*02 

8.230010+01 

9. 43717D+0I 

4.50000E-03 

3.353360*02 

8.028010+02 

6. 183810+02 

9. 144460+01 

1.048570+02 

3. 00000E-03 

3.688920*02 

8.830810+02 

6.804390+02 

1 .003890+02 

1 . 133430+02 

5 . 30000E-03 

4.024270*02 

9.633610+02 

7.422970+02 

1.097330+02 

1 . 238290+02 

6 . 00Q00E-03 

4. 339630*02 

1 .04364D+03 

8.041560*02 

1 .188780+02 

1 .36315D+02 

6 . 50000E-03 

4 . 69499D+02 

1 . 123920*03 

8.660140+02 

1.280220+02 

1.468000+02 

7 . OOOOOE-03 

5.030340*02 

1 .20420D+03 

9.278720+02 

1.371670+02 

1 .572860+02 

7.50000E-03 

5.365700*02 

1.284480+03 

9.897300*02 

1.463110*02 

1.677720+02 

8. OOOOOE-03 

3.701030*02 

1 .36476D+03 

1-031590*03 

1 .53436D*02 

1 .78238D+02 

8.50000E-03 

6.036410*02 

1 .44304D+03 

1 . 1 13450*03 

1 .646000*02 

1 . 88743D*02 

9. 00000E-03 

6.371770*02 

1 .323320+03 

1 . 173300+03 

1.737450*02 

1.992290*02 

9.50000E-03 

7  .  4  132  40 ♦  1  2 
7.413240*12 
n  .  4  1  3240* I  2 
7.413240*12 
7. 413240+12 
7.413240+12 
7  .  41 324D+12 
7.41 324D+ 1 2 
7.413240+12 
7.  41  3240+12 
7.413240+12 
7 . 41 324D+12 
7 . 4 l 324D+ 1 2 
7.413240+12 
7.413240+12 
7. 4  1  324D+12 
7.4132  4D+ 1 2 
7.41  32  4D+ l 2 
7.413240+12 


1  .  774640+13 
1 .774640+13 
1 . 77  46 4D+  l  3 
1 .774640+13 
l . 774640+13 
1.774640+13 
1 .774640+13 
1 . 77  464D+  1  3 
l .774640+13 
1 .77464D+13 
1 . 77464D+ 1 3 
t  .  774640+13 
1.774640+13 
l .774640+13 
l . 77464D+  1  3 
1.774640+13 
1 .774640+13 
1 . 77 46 4D+  l  3 
1 .774640+13 


1 . 36741D+1 3 
l . 367410* l 3 
1.367410+13 
1 . 367  4 1 D  +  1  3 
1 . 3674 1 D+  I  3 
1 . 367410+  l  3 
1.367410+13 
1 . 367  4 1 D+  1  3 
1 . 367410+13 
1 . 367410+13 
l . 367410+13 
l . 3674 1 D+  l  3 
1 .367410+13 
1 .36741D+13 
1 . 36741D+13 
1.367410+13 
l . 3674 1 D+  l  3 
1.367410+13 
1 . 367410+13 


2.021430+12 
2.02  1  43D+  1  2 
2.021430+12 
2.02143D+12 
2. 021 43D+12 
2.021430+12 
2.021 43D+ 1 2 
2.02 143D+12 
2.021430+12 
2.02143D+12 
2.02143D+12 
2.021430+12 
2. 021430+12 
2.02143D+12 
2.02143D+12 
2.021430+12 
2. 021430+12 
2.021430+12 
2.021430+12 


2 . 31793D+ 12 
2. 3 1 79 3D  +  l  2 
2.31 79  3D*  i  2 
2. 317930+12 
2 . 31793D+1 2 
2. 31793D+I2 
2.31 793D+ 1 2 
2. 31793D+12 
2. 31793D+1 2 
2.317930+12 
2 . 3 1 7930+  l  2 
2.31793D+12 
2. 317930+12 
2 . 3 1 793D+ 1 2 
2.31793D+12 
2. 31793D+12 
2 . 3 1 793D+ l 2 
2.317930+12 
2. 317930+12 


5. OQOOCE-04 
1 . 0000QE-03 
1 . 50000E-0  3 
2 . OOOOOE-Q3 
2 . 50000E-03 
3. 00000E-O3 
3 . 50000E-03 
4. 00000E-03 
4.50000E-03 
5 . 00000E-03 
5 . 50000E-03 
6 . 0OOOOE-O3 
6 . 50000E-03 
7 . 00000E-Q3 
7 . 50000E-03 
8 . OOOOOE-03 
8 . 50000E-03 
9. OOOOOE-03 
9 . 50000E-03 


TABLE  10 

PERTURBED  NUMBER  DENSITY,  SECTION  2-5.2 


T6 

T7 

T8 

T9 

T10 

R 

375140*10 

2 . 86405D+09 

1 . 46796D+  1  0 

1 . 89266D+09 

-1 .406640+10 

5 . Q0000E-04 

375140+10 

2 . 864050*09 

l  .  46796D+ 1 0 

1 .892660*09 

-l . 406640+10 

1 .OOOOOE-03 

375140+10 

2.86405D+09 

1 . 46796D+ 1 0 

1 .892660*09 

-1 .40664D+10 

1 .50000E-03 

375140*10 

2.864050*09 

l  . 46796D+ 10 

1 .89266D*09 

-1.406640+10 

2. 00000E-03 

375140*10 

2.864050+09 

1 . 46796D+ 1 0 

1.892660+09 

-1 . 406640+10 

2 . 50000E-03 

375140*10 

2 . 86405D+09 

1.467960+10 

1 .89266D+09 

-  1 . 40664D+ 10 

3. 00000E-03 

375140*10 

2.864050+09 

1 .46796D+10 

1.892660+09 

-1.406640+10 

3 . 50000E-03 

375140*10 

2.864050+09 

1 .46796D+10 

1.892660+09 

-1 .40664D+10 

4. 000008-03 

375140*10 

2.864050+09 

1.467960+10 

1 .892660+09 

-1.406640+10 

4.50000E-03 

375140*10 

2.864050+09 

1.467960+10 

1 .892660+09 

-1.406640+10 

5 . OOOOOE-03 

375140*10 

2.864050+09 

1.467960+10 

1 .692660+09 

-I .406640410 

5 . 50000E-03 

375140*10 

2.864050+09 

1.467960+10 

1.892660+09 

-1.406640+10 

6. OOOOOE-03 

3751 40*10 

2 . 864050+09 

1.467960+10 

1 . 89266D+09 

-1.406640+10 

6.50000E-03 

375140*10 

2.864050+09 

1.467960+10 

1.892660+09 

-1.406640+10 

7. 00000E-03 

375140+10 

2 . 86405D+09 

1.467960+10 

1 .892660+09 

-1.406640+10 

7 . 50000E-03 

375140*10 

2.864050+09 

1.467960+10 

1.892660+09 

•1.406640+10 

8 . OOOOOE-03 

375140*10 

2.864050+09 

1.467960+10 

1  . S9266D+09 

-1.406640+10 

8 . 50000E-03 

375140*10 

2.864050+09 

1.467960+10 

1.892660+09 

-1.406640+10 

9. OOOOOE-03 

37514D+10 

2 . 86405D+09 

1.467960+10 

1 .89266D+09 

-  1 . 40664D+1 0 

9 . 50000E-03 

TABLE  11 

PERTURBED  RADIAL  VELOCITY,  SECTION  2-5.2 


Ti  T2 

-1 . 381930*02  -4.879040*01 
-2.763860*02  -9.748090*01 
-4.14579D*02  -1.46221D*02 
-5.527730*02  -1.949620*02 
-6 . 909660*02  -2.437020*02 
-8.291590*02  -2.924430*02 
-9.673520*02  -3.411830*02 
-1.105550*03  -3.899240*02 
-1.243740*03  -4.386640*02 
-1.381930+03  -4 . 87404D+02 
-1.520130*03  -5.361450+02 
-1.658320+03  -5.848850+02 

-  1 . 7965 l D  +  03  -6.336260+02 

-  1 . 93470D  +  03  -6.823660*02 
-2.072900+03  -7.311070*02 
-2.21 1 09D+03  -7.798470*02 
-2 . 34928D*03  -8.285880+02 
-2.487 48D*03  -8.773280*02 
-2.625670*03  -9.260680+02 


T3  T4 

1.224000*02  8.84013D+01 
2.448000*02  1 . 76803D*02 
3.672000+02  2.652040*02 
4.896000+02  3.536050*02 
6.120000*02  4.420070*02 
7.344000+02  5.304080+02 
8.568000+02  6.18809D+02 
9 . 79200D+02  7.072110+02 
1.101600+03  7.956120+02 
1.224000+03  8.840130+02 
1.346400+03  9.724150+02 
1.468800+03  1.060B2D+03 
1.591200+03  1.149220+03 
1.71 3600+03  1.237620+03 
1.836000*03  1.326020+03 
1.958400+03  1.414420+03 
2.080800+03  1.502820+03 
2.203200+03  1.591220+03 
2.325600+03  1.67963D+03 


T5  R 

-9.375550+01  5 . OOOOOE-04 

-  1 . 875 1  10*02  1 . OOOOOE-03 

-2.81 267D+02  1 . 50000E-03 

-3.750220+02  2. 00000E-03 

-4. 687780+02-  2.50000E-03 
-5 . 62533D+02  3. 00000E-03 

-6.562890+02  3.50000E-03 

-7.500440+02  4. 00000E-03 

-8 . 43800D+02  4.50000E-03 

-9 . 3755SD+02  5. OOOOOE-03 

-1.031 3 1D+03  5 . 50000E-03 

-1.1 2507D+03  6. OOOOOE-03 

-1. 218820+03  6 . 50000E-03 

-1.31 258D+03  7. 00000E-03 

-1. 406330+03  7 . 50000E-03 

-1.500090+03  8. OOOOOE-03 

-  1 . 59384D+03  8.50000E-03 

-  1 . 68760D  +  03  9 . OOOC  0E-03 

-1.781350+03  9 . 50000E-03 


T6 

5.261280-01 
1 .05226D+00 
I  .578380*00 
2. 104510*00 
2.630640*00 
3. 156770+00 
3.682900+00 
4.209030+00 
4.735150+00 
5.261280+00 
5.787410+00 
6.313540+00 
6.839670+00 
7 . 36580D+00 
7.891920+00 
8.418050+00 
8.944180+00 
9 . 47031D+00 
9.996440+00 


TABLE  12 

PERTURBED  RADIAL  VELOCITY,  SECTION  2-5.2 


T7 

T8 

T9 

T10 

R 

-8.61 39 ID+OO 

-3.317270+00 

7 .539020+00 

5.76013D+00 

5 . OOOOOE-04 

-1.722780+01 

-6.63434D+00 

1 .S0780D+01 

1 . 152030+01 

1 .OOOOOE-03 

-2.584170+01 

-9.951820+00 

2.261710+01 

1.728040+01 

l . 30000E-03 

-3. 44556D+01 

-1.32691D+01 

3 . 015610+01 

2 . 30409D+01 

2. 00000E-03 

-4 . 306990+0 1 

-  1 .658640  +  01 

3.769510+01 

2 . 88006D+0 1 

2 . 50000E-03 

-5.168350+01 

-1 .990360*01 

4.523410+01 

3 . 45608D+0 1 

3. OOOOOE-03 

-6.029740+01 

-2 . 32209D+0 l 

5.277310+01 

4.032090+01 

3 . 50000E-03 

-6.891130+01 

-2.65382D+01 

6.031210+01 

4 . 608 10D+01 

4. 00000E-03 

-7.752520+01 

-2. 985550+01 

6.785120+01 

5. 184120+01 

4.50000E-03 

-8.61 391D+01 

-3.317270+01 

7.539020+01 

5. 760 130+01 

5 . OOOOOE-03 

-9.475300+01 

-3.649000+01 

8.292920+01 

6 . 336 1 4D+01 

5 . 50000E-03 

-1.033670+02 

-3.980730+01 

9.O4682D+0I 

6.9121 60+0 1 

6. 00000E-03 

-1 .  119810+02 

-4. 312450+01 

9.800720+01 

7.488170+01 

6 . 50000E-03 

- 1 . 20595D+02 

-4.644180+01 

1 .055460+02 

8 . 06418D+01 

7. OOOOOE-03 

-1.292090+02 

-4.975910+01 

1 . 1 30850+02 

8.640190+01 

7 . 50000E-03 

- 1 . 378230+02 

-5.307640+01 

1.206240+02 

9.216210+01 

8. OOOOOE-03 

-1.464360+02 

-5.639360+01 

1.281630+02 

9.792220+01 

8 . 50000E-03 

-1.550500+02 

-5.971090+01 

1.357020+02 

1.036820+02 

9. 00000E-03 

-1.636640+02 

-6 . 30282D+01 

1 .43241D+02 

1.094420+02 

9 . 50000E-03 

59 


TABLE  13 

PERTURBED  AXIAL  VELOCITY,  SECTION  2-5.2 

T 1  T2  T3  T  4  T5  R 

-3.879720-04  -9. 287560-04  -7.156330-04  -1.057920-04  -1.213090-04  5.00000E-04 

-1.551890-03  -3.715020-03  -2.862530-03  -4.231670-04  -4.852360-04  1.00000E-03 

-3.491750-03  -8.358800-03  -6.440700-03  -9.52126D-04  -1.091780-03  1.50000E-03 

-6.207550-03  -1.486010-02  -1.145010-02  -1.69267D-03  -1.940940-03  2.00000E-03 

-9.699290-03  -2.32189D-02  -1.789080-02  -2.64479D-03  -3.032730-03  2.50000E-03 

-1.396700-02  -3.343520-02  -2.576280-02  -3.808500-03  -4.367120-03  3.00000E-03 

-1.90 106D-02  -4 . 5509  ID-02  -3.50660D-02  -5.18380D-03  -5.944140-03  3.50000E-03 

-2.483020-02  -5.944040-02  -4.3800SD-02  -6.770670-03  -7.763780-03  4.00000E-03 


-3.142570-02 

-7 . 52293D-02 

-5.796630-02 

-8 . 569 1 3D-03 

-9.826030-03 

4 . 50000E-03 

-3.879720-02 

-9.287560-02 

-7. 156340-02 

-1 .057920-02 

-1 . 21309D-02 

5 . 00000E-03 

-4. 69446D-02 

-1 . 1 2380D-0 1 

-8.659170-02 

-1 . 28008D-02 

-1 -46784D-02 

5 . 50000E-03 

-5.586790-02 

-  1 . 3374  ID-0 1 

-1.030510-01 

-1 .523400-02 

-1.746850-02 

6 . 00000E-03 

-6.556720-02 

-1  -56960D-01 

-1.209420-01 

-1.787880-02 

-2.050120-02 

6 . 5Q0OOE-O3 

-7.604250-02 

-1.820360-01 

-1.402640-01 

-2.073520-02 

-2.37766D-02 

7 . OOOOOE-03 

-8.729370-02 

-2.089700-01 

-1.610180-01 

-2- 38031D-02 

-2.729450-02 

7 . 50000E-03 

-9 . 93208D-02 

-2 . 37762D-0  t 

-1.832020-01 

-2 . 70827D-02 

-3. 10551D-02 

8. 00000E-03 

-1.121 24D-0 1 

-2.684110-01 

-2.068180-01 

-3.057380-02 

-3 . 50583D-02 

8 . 50000E-03 

-  1 . 25703D-0 1 

-3.009170-01 

-2.318650-01 

-3. 42765D-02 

-3.930410-02 

9. 00000E-03 

-  1 . 400580-01 

-3. 352810-01 

-2.583440-01 

-3.81 908D-02 

-4. 37926D-02 

9 . 50000E-03 

TABLE  14 


PERTURBED  AXIAL  VELOCITY,  SECTION  2-5.2 


T6 

7 . 196770-07 
2.87*710-06 
6.677100-06 
1  .  15148D-05 
1.799190-05 
2.590840-05 
3.526420-05 
4. 605940-05 
5.829390-05 
7. 19677D-05 
8.708100-05 
1.036340-04 
1  -  2 1626D-04 
1 .41057D-04 
1.619270-04 
1.842370-04 
2.079870-04 
2.331750-04 
2.598040-04 


T7 

-1.498900-07 
-5 . 9956  ID-07 
-1.349010-06 
-2.398240-06 
-3.747260-06 
-5.396050-06 
-7.344630-06 
-9.592980-06 
-1.2141  ID-OS 
- 1 . 498900-05 
-1.813670-05 
-2. 158420-05 
-2.S3315D-05 
-2.937850-05 
-3.372530-05 
-3.837190-05 
-4. 331830-05 
-4.856450-05 
-5.41 104D-05 


T8 

-7.682560-07 
-3.073020-06 
-6.91 430D-06 
-1.229210-05 
-1.920640-05 
-2.765720-05 
-3.764450-05 
-4.916840-05 
-6.222870-05 
-7.682560-05 
-9.295900-05 
-1. 106290-04 
-1.298350-04 
-1.505780-04 
-t. 728580-04 
-1.966740-04 
-2.220260-04 
-2.489150-04 
-2.773400-04 


T9 

-9 . 905260-08 
-3 . 962 10D-07 
-8 . 91473D-07 
-1 .58484D-06 
-2.476310-06 
-3.S6S89D-06 
-4.85358D-06 
-6.339370-06 
-8.023260-06 
-9 . 90526D-06 
-1 . 19854D-05 
-1.426360-05 
-1.673990-05 
-1 .94143D-05 
-2.228680-05 
-2.535750-05 
-2.862620-05 
-3.209300-05 
-3.575800-05 


T10 

7.361630- 07 
2.944650-06 
6.625470-06 
1 . 177860-05 
1.840410-05 
2 . 6501 9D-05 
3.607200-05 
4.711450-05 
5 . 96292D-05 

7.361630- 05 
8 . 90758D-05 
1 .06008D-04 
1.244120-04 
1.442880-04 
1.656370-04 
1.884580-04 
2. 1275ID-04 
2 . 385 1 7D-04 
2.657550-04 


R 

5 . OOOOOE-04 
1 , OOOOOE-03 
I .50000E-03 
2. 00000E-03 
2.50000E-03 
3. 00000E-03 
3.50000E-03 
4. OOOOOE-03 
4 . 50000E-03 
5. 00000E-03 
5 . 50000E-03 
6. 00000E-03 
6 . 50000E-03 
7. 00000E-03 
7.50000E-03 
8. OOOOOE-03 
8. 50000E-03 
9. OOOOOE-03 
9 . 50000E-03 


3  .  35  3560+0  1  8.0280 10«-0  l 
6.7071 2D  +  0 1  I . 605600+02 
l . 006070+02  2.408400*02 
1.341420*02  3.211200*02 
1.676780*02  4.014000*02 
2.012140*02  4.816800*02 
2 . 34749D+02  5.619600*02 
2.682850*02  6.422410*02 
3.018200+02  7.225210*02 
3.353560*02  8.02801D+02 
3.688920+02  8.830810+02 
4.024270+02  9.633610+02 
4.359630+02  1.043640+03 
4. 69499D+02  1.12392D+03 
5.030340+02  1.204200+03 
5 . 36570D+02  1.28448D+03 
5.701050+02  1.364760+03 
6 . 036410+02  1 . 44504D+03 
6.371770+02  1 . 52532D+03 


6.185810*01  9.144450*00 
1.237160*02  1.828890*01 
1.855740*02  2.743340+01 
2.474320+02  3.65778D+01 
3.092910+02  4.572230*01 
3.711 49D+02  5 . 48667D+01 
4.330070+02  6.40112D+01 
4. 94865D+02  7.31556D+01 
5.567230+02  8.2300t0+01 
6.185810+02  9 . 1 4446D+0 1 
6.804390+02  1.005890+02 
7.422970+02  1.09733D+02 
8.041560+02  1.18878D+02 
8.660140+02  1.280220+02 
9.27872D+02  1.37167D+02 
9.897300+02  1.4631ID+02 
1.051590+03  1 . 55456D+02 
1.113450+03  1 . 64600D+02 
1.175300+03  1.73745D+02 


1 . 04857D+0 1  5.00000E-04 
2.097150+01  1 . OOOOOE-03 
3.145720+01  1.50000E-03 
4 . 1 9430D+0 1  2. 00000E-03 
5 . 24287D+0 l  2.50000E-03 
6.291 440+0 1  3. OOOOOE-03 
7 . 340020+0 1  3 . 50000E-03 
8.388590+01  4. 00000E-03 
9 . 437 17D+0 1  4.50000E-03 
1 . 04857D+02  5 . 00000E-03 
1.15343D+02  5.50000E-03 
1.258290+02  6. OOOOOE-03 
1 . 36315D+02  6 . 50000E-03 
1 . 468000+02  7. OOOOOE-03 
1 . 57  2860+02  7.50000E-03 
1.677720+02  8. OOOOOE-03 
1 . 78258D  +  02  8.50000E-03 
1 . 88743D+02  9. OOOOOE-03 
1 . 992290+02  9.50000E-03 


TABLE  16 

PERTURBED  RADIAL  ELECTRIC  FIELD,  SECTION  2-5.2 


T6 

■6.220770-02 
-1.264150-01 
-1.866230-01 
-2.488310-01 
-3.1 10380-01 
-3.732460-01 
-4.354540-01 
-4.976610-01 
-5.598690-01 
-6.220770-01 
-6.842840-01 
-7.464920-01 
-8.087000-01 
-8 . 70907D-01 
-9.331 I5D-01 
-9.95323D-01 
-1.057530*00 
-1.119740*00 
-1  . 181950*00 


T7 

1.295630- 02 
2.591250-02 
3.886880-02 
5.182500-02 
6.4781 3D-02 
7.773750-02 
9.069380-02 
1 .036500-01 
1 . 166060-01 

1.295630- 01 
1.425190-01 
1.554750-01 
1 .6843ID-01 
1 .813880-01 
1.943440-01 
2 . 07300D-0 1 
2.202560-01 
2.332130-01 
2.461690-01 


T8 

6.640670- 02 
1.328130-01 
1.992200-01 
2.656270-01 
3 . 320340-0 1 
3.984400-01 
4.648470-01 
5.31 254D-01 
5.976600-01 

6.640670- 01 
7.304740-01 
7.968810-01 
8.632870-01 
9.296940-01 
9.961010-01 
1.062510*00 
1 .128910*00 
1 . 195320*00 
1.261730*00 


T9 

8 . 56 193D-03 
1.712390-02 
2 . 568S8D-02 
3 . 42477D-02 
4.280970-02 
5.137160-02 
5.993350-02 
6.849S5D-02 
7 . 70574D-02 
8.561930-02 
9.418130-02 
1. 027430-01 
1 . 1 1305D-01 
1 . 198670-01 
1.284290-01 
1.369910-01 
1.4555 30-01 
1.541150-01 
I .626770-01 


T10 

-6.363270-02 
-1 .2726SD-01 
-1.908980-01 
-2. 545310-01 
-3.181630-01 
-3.81796D-01 
-4. 454290-01 
-5.090610-01 
-5 . 72694D-01 
-6 . 36327D-01 
-6.999590-01 
-7.635920-01 
-8 . 27225D-01 
-8 . 90B57D-01 
-9.544900-01 
-1.018120*00 
-1.081 76D*00 
-1 .145390*00 
-1 .20902D+00 


R 

5 . 00000E-04 
1 .000000-03 
1 .50000E-03 
2 . OOOOOE-03 
2 . 50000E-03 
3. OOOOOE-03 
3 . 50000E-03 
4.00000E-03 
4.50000E-03 
5 . OOOOOE-03 
5 . 50000E-03 
6. 00000E-03 
6 . 50000E-03 
7. OOOOOE-03 
7.50000E-03 
8. 00000E-03 
8 . 50000E-03 
9. OOOOOE-03 
9.50000E-03 


a  t 


TABLE  17 

PERTURBED  NUMBER  DENSITY,  CHAPTER  3 


-V* 


Tl 

T2 

T3 

T  4 

T5 

R 

1 

r 

7 . 66858D+ 1 2 

1 .782200*13 

1  .  34434D+- 1  3 

t  .  871  17D+12 

2.500120+12 

5.00000E-04 

7.668580*12 

1 .782200*13 

1 • 344340* 1 3 

1.871 17D+12 

2.500120+12 

1 . 000002-03 

7.668580*12 

1 . 782200+13 

1  *  34434D+  1  3 

1.871  17D+12 

2.500120+12 

1 . 50000E-03 

- 

7.66858DM2 

1.782200*13 

1 . 344340+13 

t  .87  l  17D+12 

2 . 500 1 20+ 1 2 

2 . 00000E-03 

7 . 66838D+ l 2 

1 .782200+13 

1  *  34434D+  1  3 

1 .871170+12 

2.500120+12 

2 . 50000E-03 

1 

7 . 66838D+ 1 2 

1.782200+13 

l  ■  3443404-13 

1.871170+12 

2.500120+12 

3.00000E-03 

7 . 66858D* 1 2 

1.782200+13 

1.344340+13 

1.871170+12 

2. 500120+12 

3.50000E-03 

7 . 66836D+ l 2 

1.782200+13 

t . 344340+13 

1.871 17Q+12 

2.S0012D+12 

4.00000E-03 

,  * 

7 . 668580* 1 2 

1 . 78220D+ 1 3 

1.344340+13 

1 .871 17D+12 

2.500120*12 

4.50000E-03 

L 

r 

7.668580*12 

1.782200+13 

l . 344340+13 

1  .  871170+12 

2.500120+12 

5 . OOOOOE-03 

P 

7 . 668380+ 1 2 

1 .78220D+13 

1.344340+13 

1.871170+12 

2 . 500 1 2D+ 1 2 

5 . 50000E-03 

r 

7.668380*12 

1 .782200+13 

1.344340+13 

1.871170+12 

2.500120+12 

6. 0O000E-03 

k‘- 

7.668380*12 

1 . 78220D+ 1 3 

1.344340+13 

1.871170+12 

2.500120+12 

6.50000E-03 

7 . 66838D* 1 2 

1  .782200+13 

1  . 34434D+ 1 3 

1.871170+12 

2.500120+12 

7 .00000E-03 

a 

| 

7 . 668380* 1 2 

1 .78220D+13 

1.344340+13 

1.871170+12 

2. 500120+12 

7 . 50000E-03 

7.668380*12 

1.782200+13 

1  . 34434D+ l 3 

1.871170+12 

2.500120+12 

8. 00000E-03 

*. 

7 . 66838D* 1 2 

1.782200+13 

1.344340+13 

1.871170+12 

2 . 5001 20+ 1 2 

8 . 50000E-03 

7 . 66838D* 1 2 

1.782200+13 

t .344340+13 

1.871170+12 

2.50012D+12 

9 . 00000E-03 

L 

7.668380*12 

1.782200+13 

1.344340+13 

I .871170+12 

2.50012D+12 

9.50000E-03 

/  .k,  ** 
v*\\v 


-i 


TABLE  lb 

PERTURBED  RADIAL  VELOCITY,  CHAPTER  3 


T!  T2  T3  74  T5  R 


1 . 389460+02 
■2.778930+02 
4. 168390+02 
5.557860+02 
6.947320+02 
8.336790+02 
9.726250+02 
1.11 1570+03 
1.250520+03 
1.389460+03 
1.528410+03 
1.667360+03 
1.806300+03 
1 .94525D+03 
2.084200+03 
2.223140+03 
2.362090+03 
2 . 50 104D+03 


-4.502210+01 
-9.004420+01 
-  1 . 350660+02 
-1 .80088D+02 
-2.251100+02 
-2.701330+02 
-3. 151550+02 
-3.601770+02 
-4.051990+02 
-4 . 502210+02 
•4.952430+02 
-5 . 40265D+02 
-5.852870+02 
-6.303090+02 
-6.753320+02 
-7 . 20354D+02 
-7.653760+02 
-8. 103980*02 


1 .243580+02 
2.487160+02 
3.730740+02 
4.974330+02 
6.217910+02 
7.461 490+02 
8.705070+02 
9 . 94865D+02 
1 . 1 19220+03 
1.243580+03 
1.367940+03 
1.492300+03 
1.616660+03 
1 .74101D+03 
1.865370+03 
1.989730+03 
2.114090+03 
2.238450+03 


8.531720+01 
1 . 706340  +  02 
2 . 559520+02 
3.412690+02 
4 . 26586D+02 
5.119030+02 
5.972210+02 
6 . 825380+02 
7.678S5D+02 
8.53172D+02 
9. 38490D+02 
1.023810+03 
1 . 109120+03 
1 . 194440*03 
1.279760*03 
1.365080+03 
1.450390+03 
1 .53571D+03 


-9.671320+01 
-1 . 93426D+02 
-2.901400+02 
-3.868530+02 
-4.835660+02 
-5 . 80279D+02 
-6 . 76992D+02 
-7 . 737060+02 
-8.704190+02 
-9.671320+02 
- 1 . 06  3850+03 
-1 . 160560+03 
-1 .25727D+03 
-1 . 353980+03 
-1.450700+03 
-1 .54741D+03 
-1.644120+03 
-1.740840+03 


5 . OOOOOE-04 
1 . OOOOOE-03 
1 . 50000E-03 
2. OOOOOE-03 
2 . 50000E-03 
3 . OOOOOE-03 
3 . 50000E-03 
4 . OOOOOE-03 
4 . 50000E-03 
5. OOOOOE-03 
5 . 50000E-03 
6. OOOOOE-03 
6 . 50000E-03 
7. OOOOOE-03 
7 . 50000E-03 
8. OOOOOE-03 
8 . 50000E-03 
9. 00000E-03 


2.639980+03  -8.554200+02  2.362800+03  1.621030+03  -1.83755D+03  9 . 50000E-03 


4.01 3350-04 
1 .605340-03 
3. 612020-03 
6.421360-03 
1.003340-02 
1 .444810-02 
1.966540-02 
2.568550-02 
3 . 23082D-02 
4.013330-02 
4.836160-02 
3.779230-02 
6.782570-02 
7.866170-02 
9 . 03003D-02 
1 .027420-01 
1 . 159860-01 
1 . 300330-01 
1 . 44882D-0 1 


-9 . 327 1 30-04 
-3.730850-03 
-8.394420-03 

-  1 . 492340-02 
-2.331780-02 
-3.337770-02 
-4.370300-02 
-3.969370-02 
-7.334980-02 
-9.327140-02 
-1 . 12838D-01 
-l . 3431 10-01 

-  1  -37629D-01 
-1.8281 2D-0 1 
-2.098610-01 
-2.387750-01 
-2 . 69534D-0 1 
-3.021 99D-0 1 
-3. 36710D-01 


-7.035580-04 
-2.814230-03 
-6.332020-03 
-1 . 1 2569D-02 
-1.758890-02 
-2.532810-02 
-3.447430-02 
-4.302770-02 
-5.698820-02 
-7.033380-02 
-8.313030-02 
-1 .013120-01 
-1 . 189010-01 
-  1 . 37897D-0 1 
-1 .38301D-01 
-1.801110-01 
-2 . 033280-01 
-2.279530-01 
-2.339840-01 


-9.792760-03 
-3.917100-04 
-8.81 3480-04 
-1 .366840-03 
-2.4481 9D-03 
-3.323390-03 
-4 . 79845D-03 
-6.267370-03 
-7 . 932 1 3D-03 
-9 . 792760-03 

-  1 . 184920-02 

-  1 .410160-02 
-1.634980-02 
-1 .9I938D-02 
-2.203370-02 
-2.506930-02 
-2.8301  ID-02 
-3. 172850-02 
-3.3331 9D-02 


-1 . 30844D-04 
-5.233750-04 
-1 . 17759D-03 
-2.093300-03 
-3. 271 10D-03 
-4.710380-03 
-6.41 135D-03 
-8 . 3740ID-03 
-1 . 05984D-02 
-  1 . 308440-02 
-1.5832  ID-02 
-1 .88413D-02 
-2.211260-02 
-2 . 364340-02 
-2.943990-02 
-3.349600-02 
-3.781 39D-02 
-4.239340-02 
-4.723460-02 


5 . 00000E-04 
l . 0Q000E-03 
1 . 50000E-03 
2 . OOOOOE-03 
2 . 50000E-03 
3. 00000E-03 
3 . 50000E-03 
4. OOOOOE-03 
4 . 50000E-03 
3. OOOOOE-03 
5 . 50000E-03 
6. 00000E-03 
6 . 50000E-03 
7. OOOOOE-Q3 
7 . 50000E-03 
8. OOOOOE-03 
B.50000E-03 
9. 00000E-03 
9 . 50000E-03 


3.469070*01  8.062210*01 
6.938150*01  1.612440*02 
1.040720*02  2.418660*02 
1.387630*02  3 ■ 224880*02 
1.734540*02  4.031110*02 
2.081440*02  4.837330*02 
2.428350*02  5.643550*02 
2.775260+02  6.449770*02 
3.122170*02  7.255990*02 
3.469070+02  8.06221D+02 
3.815980+02  8.868430+02 
4.162890+02  9.67466D+02 
4 . 509800+02  1 . 048090+03 
4.856700+02  1.128710+03 
5.203610+02  1 . 209330+03 
5 . 55052D+02  1 . 28995D+03 
5.897420+02  1.370580+03 
6.244330+02  1.45120D+03 
6.591240+02  1.531820+03 


6.081430+01  8.464690*00 
1.216290+02  1 . 692940+01 
1.824430+02  2.539410*01 
2.432570+02  3.38588D+01 
3.040720+02  4.232340+01 
3.648860+02  5.078810+01 
4.257000+02  5.92528D+01 
4.865140+02  6 . 77 I75D+01 
5 . 47329D+02  7.61822D+01 
6.081 43D+02  8.46469D+01 
6.689570+02  9.31116D+01 
7.297720+02  1.01576D+02 
7.90586D+02  1.100410+02 
8.51 400D+02  1.185060+02 
9.122150+02  1.269700+02 
9.730290+02  1.354350+02 
1.033840+03  1.439000+02 
1 . 09466D+03  1.523640+02 
1.155470+03  1.608290+02 


1.130990+01  5 . 00000E-04 
2.261980+01  1 . OOOOOE-03 
3.392970+01  1.50000E-03 
4 . 52397D+0 1  2. OOOOOE-03 
5.654960+01  2.50000E-03 
6.785950+01  3. OOOOOE-03 
7 . 9 1 694D+01  3.S0000E-03 
9 . 04793D+0 l  4. 00000E-03 
1.017890+02  4.50000E-03 
1 . 1 3099D+02  5. OOOOOE-03 
1.244090+02  5 . 50000E-0  3 
I . 357 1 9D+02  6. 00000E-03 
1.470290+02  6.50000E-03 
1 . 58339D+02  7. 00000E-03 
l . 69649D+02  7.50000E-03 
1 . 809590+02  8. OOOOOE-03 
1 . 92269D+02  8.50000E-03 
2 . 03578D+02  9. 00000E-03 
2 . I 4888D+02  9 . 50000E-03 


139276823847E-10 
7 08 28191391 4E-09 
2 46 6864977 93E- 07 
738177854306E-06 
2 1 9407939 393E-04 
63206852 1320E* 03 
1 9 37848 8567 0E- 01 
573899063330E*00 
171 148461 225E+02 
3086 2708042 1E>03 


.2I2743966624E-06 
. 948153203973E-06 
. 22088601 5309E- 03 
.3994858 1 32 13E-03 
.6306 1 2931330E-03 
. 914208323973E-03 
. 124832384989E-04 
. 13813891 06 40E-04 
. 368940692 1 97E-03 
- . 482883 1 82953E~03 


•329232209408E-10 
. 144833322286E-08 
. 30108386 3021E-07 
. 149851 7893 34E- 03 
. 443384221883E-04 
. 132365445078E-02 
. 393369999991E-01 
.  1  1690335 1683E+0 1 
. 34741961 1898E+02 
. 103247801 128E+04 


.  360000000000E- 1 1 
. 760000000000E- 1 1 
. 1 16000000000E-10 
. I36000000000E-10 
. 1 96000000000E- 1 0 
. 2  36000000000E- 1 0 
. 276000000000E* 1 0 
.316000000000E-10 
.336000000000E-10 
. 396000000000E-10 


TABLE  22 

PERTURBED  QUANTITIES,  CHAPTER  4 
Ez(r,t)  BQ(r,t)  Bz(r,t)  TIME 


143376253407E-06 
631619296326E-05 
2185206331 33E-03 
65349 49 46509E-02 
194229471 106E+00 
577238014437E+0I 
171546371 1 85E+03 
50980959381 1E+04 
15 I507673948E+06 
45025766 1 152C+07 


-.281 178995039E- 15 
143S98567192E-13 
-.314546648399E-12 
154121530392E-10 
-.4581212 38 718E-09 
-. 13613157 1017E-07 
-. 40462 19043 2 8E- 06 
-. 12024744731  IE-04 
-.357357162373E-03 
-. 106201089316E-01 


-. 128953332382E-18 
-. 66773908 39 49E-17 
-.2359795939 46E- 15 
-.706826801289E-14 
21010196884 4E-12 
-. 6  2441 359872  IE- 1  I 
1 85566 290305E-09 
-. 55147 46598 36 E- 08 
I63889898677E-06 
-  .  4870557 4 1426E-05 


. 360000000000E- 1 1 
■ 760000000000E- 1 1 

•  1 1 6000000000E- 1 0 

•  1 36000000000E- 1 0 

■  1 96000000000E-  1 0 
. 236000000000E- 10 
. 276000000000E- 1 0 

■  3 1 6000000000E- 1 0 
. 3360000000 00E- 1 0 
. 396000000000E- 10 


/  / 


PERTURBED  RADIAL  ELECTRIC  FIELD, Kg-M/SEC  -COULOMB 


PERTURBED  AXIAL  VELOCITY,  M/SEC 


PERTURBED  RADIAL  ELECTRIC  FIELD,  Kg-M/SF.C  -COULOMB 
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